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Abstract. Transport properties of a thermal medium determine how its conserved charge densities (for 
instance the electric charge, energy or momentum) evolve as a function of time and eventually relax back 
to their equilibrium values. Here the transport properties of the quark-gluon plasma are reviewed from 
a theoretical perspective. The latter play a key role in the description of heavy-ion collisions, and are an 
important ingredient in constraining particle production processes in the early universe. We place partic- 
ular emphasis on lattice QCD calculations of conserved current correlators. These Euclidean correlators 
are related by an integral transform to spectral functions, whose small-frequency form determines the 
transport properties via Kubo formulae. The universal hydrodynamic predictions for the small-frequency 
pole structure of spectral functions are summarized. The viability of a quasiparticle description implies the 
presence of additional characteristic features in the spectral functions. These features are in stark contrast 
with the functional form that is found in strongly coupled plasmas via the gauge/gravity duality. A central 
goal is therefore to determine which of these dynamical regimes the quark-gluon plasma is qualitatively 
closer to as a function of temperature. We review the analysis of lattice correlators in relation to transport 
properties, and tentatively estimate what computational effort is required to make decisive progress in this 
field. 
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1 Introduction 

The study of a bulk state of matter such as water be- 
gins with its equilibrium properties, such as the equation 
of state and the pair correlation function. However both 
to understand the nature of the state more deeply and 
to interpret the outcome of experiments, which necessar- 
ily involve time evolution at some level, it is essential to 
study the dynamic properties of the medium. A natural 
starting point is to focus on the response of the medium to 
long-wavelength and slow-frequency perturbations in en- 
ergy density, momentum density and conserved charges. 
This is the realm of hydrodynamics. In addition one will 
at first restrict oneself to small-amplitude perturbations, 



and expand in them to first or perhaps second order. When 
studying a medium in which quantum and relativistic ef- 
fects play a dominant role, the same approach proves use- 
ful. In this review wc focus on the state of matter ob- 
tained when heating up ordinary nuclear matter to tem- 
peratures about ISO'OOO times the core temperature of the 
Sun. If the temperature is raised gradually, the system is 
thought to consist at first of a gas of hadrons (nucleons, 
pious, kaons, . . . ). Around a temperature T — Tc wlOO- 
200MeV the hadrons 'melt', leading the system to go over 
to a phase where individual color charges can no longer be 
assigned to a unique hadron. At very high temperatures 
T ^ Tc, the color charge is transported by weakly inter- 
acting quasiparticles carrying either quark or gluon quan- 
tum numbers. This ultimate high-temperature behavior is 
a firm prediction of Quantum Chromodynamics (QCD), 
the fundamental theory that quantitatively describes the 
interactions of quarks and gluons, due to 'asymptotic free- 
dom'. The latter property of the theory implies that the 
amplitude for large-angle scattering between two energetic 
particles is small. The phase above Tc is called the quark- 
gluon plasma. It equilibrium properties have been studied 
extensively using Monte-Carlo simulations [T]. 

Two research fields provide the main motivation to 
study the dynamic properties of the quark-gluon plasma: 
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early-universe cosmology and heavy-ion collisions. In the 
latter context, a bulk state of hadronic matter is formed 
for a time of the order of 5fm/c. Therefore dynamic prop- 
erties of this hot phase are obviously important. An a 
priori open question is whether the relaxation times of 
the medium are short enough for it to be treated as being 
locally in equilibrium, thus enabling a hydrodynamic de- 
scription of the system's evolution. In the context of the 
early universe, the strongly interacting particles were in 
the quark-gluon plasma phase up to a time of the order of 
a few microseconds, at which point the transition to the 
hadronic phase took place. While the universe's expan- 
sion was slow enough to allow the quark-gluon plasma to 
equilibrate, the production of weakly interacting particles 
depends on certain dynamic structure functions, as we will 
illustrate below, with potentially observable consequences. 

Heavy ion collisions have been the experimental ap- 
proach to study the quark-gluon plasma since the mid 
1970's. In the standard interpretation of these reactions, 
a chunk is knocked out of each nucleus, and these two 
fragments undergo strong interactions and form a 'fire- 
ball', out of which many particles are produced. While 
the 'spectator' nucleons fiy down the pipeline, the de- 
cay products of the fireball are recorded in particle detec- 
tors. The typically O(IO^) particles detected in the 'cen- 
tral rapidity region' are produced with a certain distri- 
bution in the azimuthal angle (f) around the beam axis. 
The distribution is parametrized by Fourier coefficients, 

dN/d(j) (x[l + 2 J2n «n C0s{2(f>)]. 

Undoubtedly one of the most striking observations from 
the Relativistic Heavy Ion Collider (RHIC) of Brookhaven 
National Laboratory is the very large elliptic flow V2 [HISl 
ElIT]. The largest nuclei collided are gold nuclei, at a center- 
of-mass energy per nucleon of \/s/A — 200GeV. Elliptical 
flow in peripheral heavy-ion collisions is a response of the 
system to the initial, anisotropic 'almond-shaped' region 
where the interactions take place, see Fig. ([T]). The mea- 
sured size of elliptical flow suggests that the shear viscos- 
ity to entropy density ratio rj/s oi the formed system is 
smaller than ever observed in other systems . In kinetic 
theory, the shear viscosity is proportional to the mean free 
path of the quasiparticles; a short mean free path corre- 
sponds to strongly interacting matter, which is the reason 
why the hot quark matter produced at RHIC has been 
called a strongly-coupled quark gluon plasma (sQGP). 

The Pb-Pb collisions at the Large Hadron Collider 
(LHC) at CERN have very recently started to provide 
clues as to how sensitive the remarkable properties of the 
quark-gluon plasma are with respect to a change in tem- 
perature [3]. Inspite of an order-of-magnitude increase in 
the center-of-mass energy, \fslA = 2.76TeV, the charged- 
particle differential elliptic flow V2{pt) is, within the uncer- 
tainties, identical to the flow observed in the Au-Au colli- 
sions at RHIC (see Fig. [5]). Whatever the outcome of the 
ongoing analyses, the phenomenology of these collisions 
has generated a strong interest for transport properties in 
strongly coupled quantum systems f8 . 

Bounds on the shear viscosity can be extracted phe- 
nomenologically by comparing experimental elliptic-fiow 
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Fig. 1. The response of nuclear matter to an anisotropic ini- 
tial condition in configuration space, expressed as the elliptic 



flow V2{pt) per unit excentricity ehydro ~ 



where aver- 



ages are taken with respect to the number of participants in 
the transverse plane. The RHIC measurements (here by the 
STAR collaboration) are compared to the ideal hydrodynamic 
predictions. Figure from [2], presented in the review [3]. 



data to the dependence of V2 on this parameter [TUIITTIIT^ , 
as illustrated in Fig. ([3|). Viscous hydrodynamic calcula- 
tions are technically rather involved, in particular because 
a naive implementation of the leading dissipative terms 
leads to unstable solutions. However, several groups have 
come up with numerical schemes (involving second-order 
terms) that cure this problem, and the results are in good 
agreement (see the review [3].) Since hydrodynamics is a 
good description at best for a part of the time evolution of 
the fireball, there are many sources of systematic uncer- 
tainty in constraining the values of rj/ s that are consistent 
with the data. The largest seems to be the sensitivity to 
the initial conditions. Nonetheless, the confidence in the 
hydrodynamic description can be increased by confronting 
the predictions for many correlation observables jl3) to the 
experimental data. 

A conservative bound (see e.g. [3]) 
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has been inferred from the RHIC data as characterizing 
the produced medium. This number can be compared to 
the perturbative prediction, which in leading-logarithmic 
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Fig. 2. Comparison of the elliptic flow measured at RHIC 
(Au-Au collisions at -/si A = 200GeV) and at the LHC (Pb- 
Pb collisions at y/s/A = 2.76TeV) ^. 



approximation reads 



leading log g'^log{l/g) ' 



(2) 



c being a constant (see [M] and Refs. therein) and g the 
strong couphng constant. Numericahy, ^ « 2.0 for as = 
0.15 at full leading order [121 . This value of the coupling 
is relevant to temperatures well above those reached in 
heavy-ion collisions, and it is uncertain how the numerical 
value extrapolates to lower temperatures. Since the phe- 
nomenological bound is small, it has become customary 
to also compare it to the value oi rj/s in the A/" = 4 super- 
Yang-Mills (SYM) theory in the limit of infinite gauge 
coupling ^16) . 

'"^ - ' (3) 



A^=4,A=oo 
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Due to the super-conformality of the A/" = 4 SYM the- 
ory, the coupling constant does not evolve with the energy 
scale, instead it is a truly constant parameter of the the- 
ory, unlike in QCD. In fact the result is much more 
general: it holds for a large class of field theories admitting 
a dual gravity description |17j. 

There are several other phenomenological indications 
that the medium created at RHIC is strongly interacting, 
see for instance [19]. The elliptic flow of charm quarks 
through the plasma is estimated experimentally via lep- 
tonic decays of D mesons [^l?!] . The charm quarks ex- 
hibit a significant amount of elliptic flow, from which one 
infers that they kinetically thermalize on a time-scale which 
is small compared to the naively expected one, Af/T^ [22 . 
This aspect of relaxation in the quark-gluon plasma is re- 
viewed in section (|3.3p . 

Another indication comes from the study of jets, i.e. the 
passage of a fast parton through the medium which then 
fragments successively into the hadrons that are detected 
[231124] . It is found that when triggering on a jet, the total 
transverse momentum is balanced not by a second back- 
to-back jet, but is balanced instead by a large number 
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Fig. 3. The elliptic flow predicted by viscous hydrodynamics, 
for different values of the shear viscosity per unit entropy 77/5 
and for Glauber (top) and color-glass condensate (bottom) ini- 
tial conditions, compared to the measurements by the STAR 
collaboration at RHIC [THl. 



of particles of much lower momentum. This observation 
leads to the statement that the system is very 'opaque' 
to colored probes. Quantitatively it is expressed by saying 
that the medium admits a large jet quenching parameter 
q. While the latter parameter is also a transport coefficient 
of the system, no method is known at present to connect 
it to Euclidean correlation functions, and we will therefore 
not discuss it further. We refer the reader to [33] for an 
introduction to the subject. 

Before we discuss the calculation of equilibrium and 
transport properties of the quark-gluon plasma from first 
principles, we briefly remark on the relevance of transport 
properties of the quark-gluon plasma to early- universe cos- 
mology. We recall that the precise understanding of pri- 
mordial nucleosynthesis has repeatedly allowed cosmolo- 
gists to put severe constraints on many new physics sce- 
narios [261 . Controlling quantitatively the quark-gluon plasma 
era of the universe could potentially have analogous ben- 
efits. 

The QCD phase transition, being a crossover at zero 
net baryon density [27], is not expected to lead to an 
observable signal in the gravitational wave background. 
Only an effective large baryon chemical potential could 
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Fig. 4. Concentration of sterile neutrinos produced in the early 
universe [52]. The dominant uncertainty comes from the lim- 
ited knowledge of the spectral function. 



potentially turn it into a first order transition, where the 
inhomogeneous nucleation processes would make an im- 
print on the gravitational field (see e.g. [5S] and Refs. 
therein) . The probability per unit volume per unit time to 
nucleate a low-temperature phase bubble out of the high- 
temperature phase is proportional to a linear combination 
of the shear and bulk viscosities in the high-temperature 
phase . In this way the transport properties would have 
an impact on cosmological observables. 

A sterile, right-handed neutrino of mass m « lOkeV 
has been proposed [30] as a candidate for warm dark mat- 
ter in the context of the :^MSM. In a given scenario the 
production of these particles happens to peak around tem- 
peratures of 0.1GeV..lGeV. Their distribution is then de- 
termined by the spectral functions of active neutrinos in 
a thermal QCD bath. The latter in turn can be expressed 
in terms of the thermal vector and axial-vector spectral 
functions of QCD [31]. Given the range of temperatures 
in which these spectral functions need to be calculated, 
non-perturbative methods are required for a precise and 
reliable result, see Fig. (g]) reproduced from [35] ■ The ex- 
ample is specific, but it is likely that other particle physics 
models with cosmological implications can be constrained 
if the transport properties of QCD are determined reliably. 

In order to understand the transport properties of the 
quark-gluon plasma, some prior familiarity with its equi- 
librium properties is required. We now summarize some 
of the essentials. The most recent published lattice QCD 
result for the entropy density as a function of temperature 
is depicted in Fig. ^ (^33 ; see also [IJ). The rapid rise 
between ISOMeV and 250McV in the number of degrees 
of freedom exerting pressure is clearly visible. The ratio 
s/T^ then saturates and only very slowly tends towards 
the asymptotic 'Stefan-Boltzmann' value of 
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One of the defining properties of a plasma is the existence 
of a finite correlation length over which (chromo-)electric 
fields are screened. Its inverse defines the Debye screening 
mass. In a weakly coupled plasma of quarks and gluons, 
the corresponding mass is calculable in perturbation the- 

rnl^g'{f + f)T\ (5) 

More generally, a widely accepted definition of m/^ is that 
it is the smallest mass associated with a screening state 
odd under Euclidean time reversal [51] , which corresponds 
to the product of time reversal and charge conjugation in 
Minkowski space. A unique feature of non-Abelian plas- 
mas is that magnetic fields are also screened (unlike in the 
Sun), and at high temperature the corresponding light- 
est screening masses are parametrically of order g^T. The 
smallest screening masses of states with a non-trivial fla- 
vor quantum number are close to 27rT, which can be un- 
derstood in terms of the minimal frequency ttT associated 
with a fermion field in the imaginary-time formalism. Fi- 
nally, as Eq. ([2]) illustrates, the mean path of a quark or 
gluon quasiparticle before it changes its direction of flight 
by an 0(1) angle is {g'^T)~^ [33]. These length scales asso- 
ciated with different physical phenomena are parametri- 
cally different at sufficiently high temperatures, but there 
are indications that at temperatures that can be reached 
in heavy-ion collisions, these scales are overlapping and 
even inverted in some cases, see the discussions in [351137] . 

Without the RHIC and the LHC heavy-ion experi- 
ments, it would currently be unclear which of the strong or 
weak coupling paradigms is more appropriate to describe 
the quark-gluon plasma at temperatures of a few Tc- The 
answer must be expected to depend to some extent on the 
physical quantity. The convergence of perturbative calcu- 
lations at T ss 2 — 3Tc (see for instance [35]) is in many 
cases poor, and these calculations are in tension with phe- 
nomenological bounds such as Eq. ([Ij. It is therefore de- 
sirable to develop non-perturbative methods to compute 
transport properties of the quark-gluon plasma starting 
from the microscopic theory, QCD. This constitutes one 
of the central themes of this review. For illustration, we 
describe the recipe to calculate the shear viscosity in the 
next paragraph. The origin of this recipe will be described 
in detail in section p.2p . 

The shear viscosity t] parametrizes how efficiently the 
momentum of a layer of fluid (assume the momentum to be 
in the plane deflned by that layer) diffuses in the direction 
orthogonal to the momentum. The diffusion constant is 



■P 



(6) 



1) 



• ANrN, 



(4) 



where e -I- p is the enthalpy density. The time-evolution of 
transverse momentum is related to the transport of shear 
stress, Txy, through the momentum conservation equation, 
doT°y + d^T^y = 0. A second key idea is that the re- 
sponse of the fluid to externally applied shear stress is en- 
coded in the equilibrium ensemble, namely in the way the 
thermal fluctuations of shear-stress dissipate in real time. 
This is the outcome of the linear response formalism, re- 
viewed in section p.2p . The time-correlation is encoded in 



Harvey B. Meyer: Transport Properties of the Quark-Gluon Plasma 



5 



-1 — I — I — I — I — I — I — I — I — [— 
' SB' i 




I I I 



20 



200 



400 600 
T[MeV] 



800 



1000 



Fig. 5. The 2+1 flavor QCD entropy density in units of T'^ as 
a function of temperature |33) . 



the commutator-defined correlator {\Txy{t, x), Txy{0)]). Its 
Fourier transform is given by the spectral function p{lu, k). 
One can then derive a 'Kubo formula', which gives the 
shear viscosity in terms of the low-frequency part of the 
spectral function, 



1- ,• 

1] = TT lim lim . 



(7) 



Furthermore, the spectral function p is formally related to 
the Euclidean correlator by analytic continuation. 



np{uj) = Im 



(8) 



where a;„ — 2'KTn is one of the Matsubara frequencies. 
Alternatively, in the mixed time-|-spatial momentum rep- 
resentation, the connection is given by an integral trans- 
form, 



GE{t,k) = / dwp(w,fe) 



cosh 



sinh /3cj/2 



(9) 



where the kernel is the thermal propagator of a free scalar 
field. The computational recipe is thus formally straight- 
forward: calculate the Euclidean correlator Ge, determine 
the spectral function by analytic continuation, and finally 
read off the shear viscosity from the slope of p at the origin. 
In practice however there are very significant difficulties 
in carrying out this program. The Euclidean correlator is 
normally only obtained in some approximation, either be- 
cause it is expanded in powers of the coupling constant, or 
because it is obtained from Monte-Carlo simulations. In ei- 
ther case, an analytic continuation is not possible without 
further information or assumptions. In particular, while 
the Euclidean correlator can be expanded in powers of 
the strong coupling, the parametric dependence of the vis- 
cosity ^ makes it clear that some form of resummation 
of the perturbative series is required. In general it is not 
known how to implement the resummation in this way, 
and an effective kinetic theory treatment has been em- 
ployed instead in 14' . In the representation ©, which is 



the form most commonly employed when the Euclidean 
correlator is determined by Monte-Carlo simulations, it is 
clear that the smoothness of the kernel implies that the 
fine features of p{ijj) (particularly on scales Auj < T) are 
encoded in the correlator in a way which is numerically 
very suppressed. 

We wish to illustrate the latter point somewhat more 
concretely. In the limit of zero-temperature, the kernel in 
^ simply becomes e~"', so that the Euclidean correlator 
is simply the Laplace transform of the spectral function. 
Since Lorentz symmetry is restored in the limit T — 0, the 
problem is simplified by going to four-momentum space. 
We will consider the case of the electromagnetic current 
correlator, j^{x) = ^ujf^u — ^d'j^d — ■^sj^s + . . . Cur- 
rent conservation and Lorentz invariance implies that a 
single function of four-momentum squared determines all 
correlation functions of the electromagnetic current. 



d'a; {j^{x)j,m e'"-^ = {q^q, - q^g^,) n{q^). (10) 



In Euclidean four-momentum space, the spectral repre- 
sentation of the vacuum polarisation n{Q^) reads 



77(0) - n{q^ 



ds 



s(s + q^) 



(11) 



Via the Optical Theorem, the spectral density is accessible 
to experiments 



7rp(s) 



47ra(s 



■o'tot{e~^e — ^ hadrons) 



a(s) 



^had(s). 



(12) 

where a is the QED coupling and the first equality holds 
up to QED correctionfQ. To illustrate the different charac- 
ter of a retarded correlator in the time-like and the space- 
like region, we reproduce in Fig. IH] the quantity i?had, as 
determined by experiments spanning several decades, as 
well as the corresponding vacuum polarization n{Q^) in 
the space-like region. While many vector resonances show 
up in i?had (and therefore the spectral function), the vac- 
uum polarisation in the space-like domain is a very smooth 
function. In the absence of auxiliary information, extract- 
ing the existence, position and width of each resonance 
from n{Q^) in the space-like region is therefore a numer- 
ically ill-posed problem. 

Returning to the physics of the quark-gluon plasma, 
it is not known precisely how far up in temperatures the 
sharp features of the vacuum spectral function depicted 
in Fig. ^ survive. However, in the case of conserved cur- 
rents, the functional form of the spectral function in the 
region uj ^ Tmax: where T,„ax is the longest relaxation time 
in the system, is completely determined by the continu- 
ous global symmetries of the theorjQ. The effective theory 



Neglecting the muon mass and up to QED radiative cor- 
rections, the quantity -Rhad can be expressed as iihad(s) ~ 

a(e'^e~ — ^hadrons) i ■ i • • j_ c ■ j_ i 

^ which IS convenient from an experimental 



a (e 



point of view [39) . 

^ At a phase transition, additional slow modes associated 
with an order parameter (possibly of a discrete symmetry) may 
play an important role. 
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Fig. 6. The R ratio, Eq. (|12|l . determined experimentally 
(Fig. from |39]), and the subtracted vacuum polarization 
n{Q^) = n{0) — n{Q^) in the space-like domain calculated 
using Eq. (pTIT2)l i40]. 

that describes this kinematical regime is hydrodynamics. 
This aspect is reviewed in section p.2p . The universality 
of these features means that they can be used as an ansatz 
for the spectral function in the small-w region. 

Beyond these universal features, the stronger assump- 
tion that the conserved charges are carried by quasipar- 
ticles leads to additional distinctive features of the spec- 
tral function around uj — 0. The latter can be predicted 
in detail by using kinetic theory. However, thanks to the 
AdS/CFT correspondence, we now know of at least one 
theory where these features are absent (the strongly cou- 
pled A/" = 4 SYM theory). We invite the reader to com- 
pare figures (TTTl black dashed curve) and (jl41 red curve), 
which illustrate the difference between the functional form 
of the shear stress spectral function in a kinetic descrip- 
tion (Boltzmann equation) and in a theory which obvi- 
ously does not admit any quasiparticles. The most impor- 
tant qualitative question about the quark-gluon plasma is 
whether, at the temperatures explored in heavy-ion colli- 
sions, it contains quasiparticles as assumed in kinetic the- 
ory. The alternative possibility is that there are instead 
no poles near the real axis in the retarded Green's func- 
tion, other than those dictated by the symmetries of the 
theory; this possibility is realized in the strongly coupled 
SYM theory. To answer the question from first principles, 
we need to study the spectral functions of conserved cur- 
rents in the small-frequency regime, in other words we 
need to study the transport coefficients. 

This review is organized as follows. We start in sec- 
tion (Sec. [2]) with a review of the analytic properties of 
thermal correlators, detailing in particular the relation be- 
tween Euclidean and real-time correlators and giving the 



general tensor structure of correlators in the channels of 
interest (Sec. 12.51) . We also describe the derivation of Kubo 
formulae and of the hydrodynamic predictions for spec- 
tral functions (Sec. 13. 2p . as well as the perturbative and 
kinetic theory treatments (Sec. l3.3|) . These results are con- 
trasted with those obtained in strongly coupled theories 
via the AdS/CFT correspondence (Sec. 13. 6|) . An overview 
of the state-of-the-art in lattice calculations of Euclidean 
correlators is given in Sec. as well as of the scaling 
of their computational cost; we discuss the electric con- 
ductivity, the heavy-quark diffusion and shear & bulk vis- 
cosities. The analytic continuation problem is discussed 
from a numerical perspective in Sec. ([5]), where we show 
the advantages and disadvantages of linear and non-linear 
methods. We finish with a summary in Sec. 

There are a number of subjects closely related to those 
treated in this review, which are only briefly mentioned 
here. The extraction of a shear viscosity estimate from 
heavy-ion data is reviewed in [S]. The thermodynamic 
properties of the quark-gluon plasma [31 , 42 , 43 , 44 ,3S] have 
been reviewed recently [1] . The behavior of spatial correla- 
tion lengths is discussed in [3S1IT] . For recent developments 
in hydrodynamics we refer to }47j . For introductions to 
relativistic hydrodynamics and the gauge/gravity duality, 
see [^H^ISO] . Finally, the physics of shear viscosity and 
its perturbative calculation are succinctly reviewed in |51) . 

2 Definition and Properties of Thermal 
Correlators 

In this section we review the different correlators rele- 
vant to calculations of transport properties. We introduce 
the Minkowski-space correlators and the Euclidean-space 
correlators, and describe how they are interrelated. Their 
spectral representation is given, from which the dispersion 
relation expressing the retarded correlator in terms of the 
spectral function is given. Finally, we summarize the ten- 
sor structure of the vector current and energy-momentum 
tensor correlators. In order to facilitate the initiation of 
the reader not accustomed to finite-temperature field the- 
ory, we begin with a reminder on vacuum correlators and 
then generalize the definitions for a thermal ensemble. 

2.1 Vacuum correlation functions 

We start with vacuum correlation functions, working in 
the metric ( — h ++). The time-evolution operator 

U{t) = e-'"* (13) 

is used to define the operators in the Heisenberg picture. 

Ait) EE iUit))-'Amit). (14) 

We recall the time-ordering operation for operators, 

T{A{h)B{t2)) = A{h)B{h) 9{ti-t2)±B{h)A{h) 9{h-h), 

(15) 
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where the + (— ) sign apphes to bosonic (fermionic) oper- 
ators. Time-ordered correlation functions, 



dte"^' {0\T{A{t)B{0))\0), (16) 



have simple Lorentz invariance properties. For instance, 
the 'Feynman' propagator of a free scalar field reads 



(0|r(<^(i,a;)<^(0))|0) 



(17) 



Up to the factor i, the Euclidean space propagator has the 
same expression, but the scalar products are to be inter- 
preted with the (-I- + ++) metric, and the ie prescription is 
unnecessary, since the integration region is well separated 
from the poles Pq+ = —m^. 

The commutator [A{t), B{0)] vanishes outside the light- 
cone in a relativistic causal theory. In the free scalar field 
theory, 



The correlator 



d^p 



pipx 



(27r)^ ipo + ^e)2 _ p2 _ ^2 • 
(18) 



Crn-* / die-*(0|[A(0,S(0)]|0), (19) 
Jo 

will play an important role in the following. In the vacuum 
it is related to the time-ordered correlator which is fa- 
miliar from time-dependent perturbation theory, through 



iG^^ {cj) 



iG 



Reuj > 
Reuj < 0. 



(20) 



Time-ordered products of fields appear when express- 
ing the time-evolution operator, and in particular the 5- 
matrix, in terms of the interaction-picture Hamiltonian. 
It is then natural to work with time-ordered propagators. 
Retarded correlation functions occur typically when cou- 
pling a field to an external source, which is switched on 
only for a finite time interval. A classic example is the 
calculation of the bremsstrahlung emitted by an electron 
passing near a nucleus. In the following we will deal mainly 
with this second type of situations, where the response 
comes however from a thermal medium rather than a sin- 
gle particle. 



2.2 Finite-temperature correlation functions 

In order to describe quantum statistical properties, the 
vacuum expectation values of operators are replaced by 
averages over an ensemble of states described by a den- 
sity matrix p. The finite-temperature, equilibrium density 
matrix is 

p^^e-^". (21) 



with Z such that Tr{/5} = 1. At finite temperature, the 
Wightman correlation functions are defined as 



G^^{t)^TT{pA{t)BiO)}., 
G^^(i)^Tr{/5B(0)A(t)}, 



(22) 
(23) 



Obviously, due to time-translation invariance of the equi- 
librium density matrix. 



Gi^{t) - G^^'i^t). 
The reality property 

G^'^'(t) -G'f^(-r)* 



(24) 



(25) 



is also easily proven. Last but not least, the definition 
implies the Kubo-Martin-Schwinger (KMS) relation. 



G^^{t)=G^^{-t~t(3). 



(26) 



In fact, it can be shown that this property uniquely char- 
acterizes the equilibrium density matrix. An obvious con- 
sequence of Eq. ([25)1 is that the thermal correlation func- 
tion of an anticommutator is periodic in the imaginary 
time direction. 

The expectation value of a commutator, 

G^^it) = ^Tl■ {p[A{t),Bm} = z (G^^(t) - G^^(t)) , 

(27) 

is physically important. The commutator vanishes outside 
the light-cone, a reflection of the causality of the theory. 
Using Eq. (l24l) and ((25)) . one obtains the properties 



(28) 
(29) 



G^^(-t) = -G^^(t), 
G^'^'(t) = G^^(t*)*, 



The Euclidean correlator is deflned as 

Gi^{t)^Gi^{-it). (30) 
As a special case of I^E^, it obeys the relation 



Gi-^(/3-i) = G^^(t). (31) 

The Fourier transform of G^^ [t) deflnes the so-called 
spectral function El 



1 r+oo 



2m 



G'^^(t). 



It enjoys the symmetry properties 



p^'^\u)=p--{L.r. 
The integral tranform over the positive half-axis 



Gji^{uj) ^ / d<e'"*G^^(t) 



(32) 

(33) 
(34) 

(35) 



Z 



^ We use a small 'hat' on the density matrix to distinguish 
it from the spectral function. 
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is analytic in the half complex plane Ini(aj) > 0. It is 
called the retarded correlator and plays a central role in 
linear response theory, as we shall see shortly. The Eu- 
clidean correlator can be expressed as a Fourier series on 
the interval <t < (3, 



G 



dte"^'*GEit) ., 



(36) 



(37) 



where ujg = uim ■ ^ = 2ttT i and we have dropped the label 
specifying the operators A, B. 

For real w and in finite volume, p is a distribution, 
namely a discrete sum of delta functions. However, in the 
infinite- volume limit of an interacting theory, it is typically 
a smooth function, except possibly at a finite number of 
points. The most useful relations among the correlators 
are those that survive the infinite- volume limit. 

Using the KMS relation one easily shows 



^_e-*({A(0,S(0)}±) = (l±e-^") 



(38) 



°° di 



e*"*G^^(t), 



whence it follows that the spectral function is related to 
the anticommutator-defined correlator via 



tanh(^) 



^e-*({A(t),B(0)}+) 



(39) 



This relation is known as the fluctuation-dissipation theo- 
rem. Since the retarded propagator G^^ can be expressed 
in terms of p^^ via a dispersion relation, it so follows that 
the knowledge of the anticommutator-defined correlator 
also determines G^^. 



2.3 Spectral representation and relations between 
correlators 

Using (1211) and (UHl), one shows that 

P^^M = ^(G^^,^M-Gr^'H*). (40) 

In particular, for B = A\ the spectral function is identi- 
cally related to the imaginary part of the retarded corre- 
lator. 



p^^^(c.) = -ImGr'(^)e 

TT 



(41) 



The spectral representation of the correlators is use- 
ful to prove further relations between the Euclidean and 
real-time correlators. Inserting two complete sets of en- 
ergy eigenstates in the definitions of G{t) and Gsit), one 
obtains 



G^^it) = |^A_i3_e-^(^"+^")/2 (42) 



.sinh(e%^)e- 



Gs^it) - — '^AmnB 



(43) 



where we employ the notation 

Er,ra = E^-E^, A„™ = (n|A(t = 0)|m) . (44) 
The expression for the retarded correlator reads 



2 X > Ai^ijiB^ 



(45) 



-/3(£;„+£;„)/2 



sinh( 



and the spectral function for a; e M takes the form 
2sinh^ ~ Z Z^^™"-""™ 

2 m.,n 

■ e-'^(^"+^'")/2<5(c.-i?„„). 

For B ^ A^ , this formula can be obtained by taking the 
imaginary part of Gfi{uj + ie) and using the standard rep- 
After 



resentation of the delta function, S{uj) 
some algebraic manipulations, one finds for the Fourier 
coefficients (1571) 



sinh(;3£;„„/2) — 



Br, 



Formulae (|^5t and (|Tf)) show in particular that 



iujf, + Enm 

(47) 



GR{iuJt) = G 



^0. 



(48) 



Thus the frequency-space Euclidean correlator is the ana- 
lytic continuation of the retarded correlator; the situation 
in the complex lu plane is illustrated in Fig. ([7]). For ^ = 0, 
one has to remember that the retarded correlator was ini- 
tially defined for Im(a;) > 0. Therefore its value at the 
origin can only be defined as lime_>.o_^ Gulie). This value 

differs from G^^ if the spectral function is singular at the 
origin. Specifically, 

^ = A6{uj) + finite ^ G^"^ - lim Gn{ie) ^ A. 

UJ £-S-0+ 

(49) 

Finally, it is easy to verify, using Eq. (|43t and (|46)) . that 
the configuration-space Euclidean correlator can be ob- 
tained from p via 



Gi^{t) + Gr{l3-t) 



AB, 



.r., coshw(l-t) 
-^-^ sinh/^L/2 ' 



' — OO 

Gi^{t)-Gi^{P~t) 



(50) 



(51) 



°, ,r., sinhcj|-t 
^ sinh/3a;/2 



These formulae show that the part of the Euclidean cor- 
relator that is even (odd) around i = /3/2 is in correspon- 
dence with the odd (even) part of the spectral function. 
In particular, B = A implies that the Euclidean correlator 
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is symmetric around t ~ j3/2 and that the spectral func- 
tion is an odd function, while B = implies that both 
(t) and p'^^'' (w) are real. In most cases of interest, 
B ~ A ^ A\ so that both properties are satisfied. 

A simple calculation based on the explicit expression 
for the spectral density (j46p . or alternatively a contour 
integration, then leads to the Kramers-Kronig relation 

du piuj) ^ ^ — 2 — Gn{iuJi) , Vwi > 0. (52) 



Starting from (I52p . it is easy to derive by recursion the 
following equalities. 



J2n 



2ri+l 



(it^i) 



dujp 



(2n+l) 



(^)- 



LJl^O (2ri+l)/ 



TTP^-->{UJ^0). (53) 

Thus at the origin, the even derivatives of the retarded cor- 
relator along the imaginary axis are given by integrals over 
the spectral density, while the odd derivatives are equal 
to the corresponding derivatives of the spectral function, 
see the bottom panel of Fig. ([7]). 

2.4 Properties of the Euclidean correlator 

We have seen that in configuration space, the Euclidean 
correlator is defined as the analytic continuation of the 
Wightman correlator, while in frequency space it turns out 
to be the analytic continuation of the retarded correlator. 
Here we discuss further relations between Euclidean and 
real-time correlators. 

In our subsequent considerations, it will prove useful to 
express the Euclidean correlator in terms of the real-time 
correlator G(t), in order to analyze the contribution of 
hydrodynamic modes. The relation between the Euclidean 
correlator and the real-time correlator is 

sinh(27rTu) 



GEit) = T 
The moments 



du G(u) — , ^ , 
^ ' cosh(27rTM) 



cos(27rTi) ' 



(54) 



" Jo sinh^ 



(55) 



(-1)"/ dtGit) — t^nh{7rTt) 



are thus expressed in terms of integrals of G{t). 

The positivity of p{uj)/lo implies the positivity of cer- 
tain linear combinations of the Euclidean correlator. For 
instance. 



^(0) _ 2 



3Gi°)-44^)-f4^^) = 12c.^ 



duj p{lo) 
^ duj 

oo ^ 

p{uj) 



(56) 



(57) 



real 



O 
< 



o 



/slope= -r| / AGr(O) 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



(B/23IT 

Fig. 7. Top: Analytic structure of the retarded correlator 
Gr{u)) in the complex uj plane. At the circled values oiu on the 
positive imaginary axis, it coincides with the Euclidean corre- 
lator (Eq. I48|l . Bottom : At the origin, the odd derivatives of 
Gr along the imaginary axis coincide with the corresponding 
derivatives of ImGfi along the real axis. These odd derivatives 
are given by transport coefficients (section |3.2|I . The figure il- 
lustrates a case where one subtraction has been made to make 
the dispersion integral converge. 



We note that relations ([Ml) and ([56H57|) only hold if the 
integrals converge in the ultraviolet. In the cases of in- 
terest, at least one subtraction (for example between two 
different temperatures) is necessary. 

In a number of contexts, it is natural to study the 
difference between the finite-temperature and the zero- 
temperature spectral density. For that purpose it is conve- 
nient to define the so-called reconstructed Euclidean cor- 
relator by 



G'i-{t,T-T')= / du:p{Lo,T') 



coshw('| 



The identity (0 < t < ;3) 

coshw(^ — t) 
sinha;/3/2 



sinhaj/3/2 



(58) 



(59) 



allows us to derive in particular the exact relations 

GT{t.T;Q)^Y.GE{\t + mplT = Q), (60) 

G'^^^it, T; \T) = GE{t, \T) + Gb(/3 - t, \T) . (61) 
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2.5 Tensor Structure of Thermal Correlators 

The tensor structure of current correlators and energy- 
momentum tensor correlators was worked out systemati- 
cally in [32] ■ We simply report the results in the notation 
of section (|2.1|) . Setting 



(62) 



and defining B{t) analogously in terms of J'^ , we denote 
their retarded correlator by G'^{k), with k = (fco, k). Due 
to conservation of the current, 



Then one defines the tensors 



P„ 



fc2 



0, P^, =0 P, 



ij ~ "t-j 



(63) 

(64) 
(65) 



The current-current correlator then has the generic tensor 
structure 



Gfl,^.(fc) = Pj, 77^ (fco, fc') + iI^(fco, fe' 



(66) 



Thus there are two independent functions of two variables 
(fco)fc^) that characterize all possible current correlators 
at finite temperature. When = = U, one recovers 
the familiar vacuum form of the correlators, which are 
characterized by a single function of k^. 

For the energy-momentum tensor, one proceeds simi- 
larly. Setting 



(67) 



and defining B{t) analogously in terms of T"^, we denote 
their retarded correlator by G^'^"^(fc), with fc = (fco,fe). 
At zero temperature and in D spacetime dimensions. 



G 



R 



(fc) = pt'''P»l^GB{k^) + H'"''°''^^Gs{k^), (68) 



The tensor H is the polarization tensor of a spin-two par- 
ticle. The correlators of the energy-momentum tensor are 
thus characterized by two independent functions. At finite 
temperature, the tensor H splits into three independent 
tensors, which one may choose as follows 

(70) 

oL , ^P^^Pl^ (72) 



There are thus five independent functions that specify all 
correlators of T^^. The function Gi describes the shear 
channel, G3 describes the tensor channel, while G2, Gl 
and Ct describe the sound channel and the bulk channel. 
When Cl = Ct and Gi = G2 = G3, one recovers the 
zero-temperature form 



3 Retarded correlators: physical significance 
and analytic calculations 

In the previous section, we introduced the definitions of 
the Minkowski and Euclidean correlators, and described 
their basic analytic and tensorial structure. Here, using 
linear response theory, we want to show why the retarded 
correlator Gii{Ld,q) is physically important. We then go 
through the most important analytic techniques to calcu- 
late the thermal spectral functions in different regimes. We 
describe how hydrodynamics makes a prediction for the 
low frequency and momentum dependence of Gfl(a;,q), 
and thereby establish the Kubo formulae which allow one 
to compute transport coefficients ab initio (section 13. 2[) . 
When the hydrodynamic regime arises from the many- 
body dynamics of weakly interacting quasiparticles, ad- 
ditional structures can be predicted in the spectral func- 
tions using kinetic theory (sections 13.31 13. 4p . We also re- 
view the predictions of perturbation theory and the opera- 
tor product expansion at high frequencies in section p.Sp . 
Some important results on spectral functions obtained us- 
ing AdS/CFT results are summarized in p.6p . Thermal 
sum rules provide useful constraints on the spectral func- 
tions, and those are discussed in section p.7|) . Finally, cer- 
tain aspects of spectral functions in finite volume (where 
lattice QCD simulations take place) are pointed out in 
section 



3.1 Linear response theory 

When studying dynamical properties of a thermal medium, 
it is natural to start with small amplitude departures from 
equilibrium. The linear response of the system to 'slow' 
perturbations is sensitive to the first order transport co- 
efficients, as well shall see in the next section. 

We assume that the time-dependent perturbation of 
the system by an external classical field f{t) coupling to 
operator B is described by a Hamiltonian of the form 



Hf{t) = H-f{t)B{t), 



(73) 



The general tensor form of the correlators in a rotationally 
invariant average over states is 

GR,^.Mk) = {pl^A + -2iPl^.Ptp+p!:uPlp))CT{k^,k^) 

+ {P^^i^Pafl + h^PjivPtl}+Pfi^uPll3))CL{ko,k'^) 

+S^i,^apGi{ko, k^) + Q^i/,Q/3G2(fco, fc^) + Lf^^^afjGziko, k^). 



H being the Hamiltonian of the unperturbed system. The 
perturbation leads to a 'response' of physical quantities, 
i.e. a change in their expectation values with respect to 
the unperturbed ensemble. The evolution equation of an 
operator A is given by 



[Hfit),A{t)] 



(74) 



One then finds that to linear order in /, the expectation 
value of A in the perturbed system minus its unperturbed 



Harvey B. Meyer: Transport Properties of the Quark-Gluon Plasma 



11 



value is 



S{A{t))^{A{t))f-{Am (75) 



Equation (j75p is the master formula of linear response 
theory. It shows that the correlator G^^ determines the 
response of an observable A to a time-dependent external 
field that couples to B. 

A source term of the form 



fit) = e''9{-t)fo 



(76) 



is often adopted to study how the system relaxes back 
to equilibrium after having been perturbed adiabatically. 
The static susceptibility is defined as the expectation value 
of ^ at t = 0, 



From ([75)) . it follows that 



,AB 



dte"^'G^^(i) =G^J,^(^e) 



(77) 



(78) 



Looking back at Eq. the zero-frequency Euclidean 

correlator is thus equal to the static susceptibility, barring 
any delta function at the origin in the spectral function. 

Integrating both sides of ([75]). da;e*"*(.), one ob- 
tains for the adiabatic perturbation (j76p 



G 



i^{u;)fo = {SAiO))f+zu; 



dte''^'{6A{t))f. (79) 



This formula shows that the relaxation of observable A 
back to its equilibrium value determines the retarded cor- 
relator G^^{uj). Since the late-time relaxation is described 
by hydrodynamic evolution, this equation can be exploited 
to obtain a prediction of the small-w functional form of 
Gr, thereby establishing the famous Kubo formulae for 
the transport coefficients (section (j3.2p ). 

In the special case where i? is a conserved charge, 
i.e. [H, B] = 0, the situation described by e = simply 
corresponds to looking at the system with a small 'chemi- 
cal potential' Jq. Then S{A{0)) is given by the thermody- 
namic derivative 

S{A{Q)) = (3fo {{Bmm - {A){B)) , (80) 

which, by comparison with (j78p . shows that it is natural 
to assign to the retarded correlator at the origin the value 



Gr{0)=H{A{0)BiO))-{A){B)) 



(81) 



3.2 Hydrodynamic predictions and Kubo Formulas 

There are recent excellent reviews on the modern point of 
view on hydrodynamics [TTllHlfiB] . Therefore our goal here 
is simply to describe the important steps that lead to the 
Kubo formulae. Hydrodynamics can be thought of as a 
low-energy effective theory: 



— it makes predictions for the low-momentum behavior 
of correlation functions. 

— it has unknown 'low-energy constants', which are called 
the transport coefficients in this context: they are the 
coefficients of a derivative expansion. 

— their values can be determined by a matching proce- 
dure with the underlying quantum field theory. 

— the number of coefficients grows with the order of the 
expansion. 

For illustration we choose as a simple example the case 
of the diffusion of a massive particle species in a thermal 
medium. In the hydrodynamic treatment, the conserva- 
tion of the particle number is expressed by the classical 
equation 

dtn + W -j = 0. (82) 

Second, a phenomenological 'constitutive equation' must 
be introduced to close the system of equations for {n,j). 
To leading order in gradients, it takes the form of Fick's 
law, 

j = -DVn. (83) 

The conjugate variable to particle density is the chemical 
potential. The particle number can be perturbed by 



Hf^ = H — J dx /i(t, x)n{t, x), 
fi(t,x) = ^(a;)e^*6'(-t), 



(84) 



(see Eq. ^ with f ^ n and B ^ n). 
We define 



ri(a;, k) — i e 
'o 



M(fe) 



dx e 



n{t, x) 

-f4x). 



(85) 
(86) 



The conservation equation ([5^ becomes, with the help of 
Fick's law, the diffusion equation 

dtn{x) = DV^n{x). (87) 

Upon Fourier transformation, the solution reads 

n{0,k) 



n{u}, k) = 



-ico + Dk 



The initial condition n(0, k) is determined by the static 
susceptibility (see Eq. (|77l78p ). 

n{0,x) = /o°°die-'7da;'G""(t,a;-a;')M(a;'), (89) 



G""(t,a;) = i([n(t,a;),n(0)]), 
or equivalently. 



n(0,fc) =Xs(fc)Mfc), 
X.(fe)=/„°°dte-^*G""(i,fe) 

G™(i, fe) = / da;e-*''' '^G""(i, x) 



(90) 



(91) 
(92) 

(93) 



This completes the prediction for the long-wavelength, 
late-time evolution of the particle number density. In Fourier 
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space, the relaxation rate of a particle density pertur- 
bation of wave-vector k is encoded in the position of a 
pole in the lower-half of the complex uj plane, namely 



Wpolc(fc ) 



'iDk . The residue of the pole is given by 



the initial condition of the relaxation process, which is 
determined by the static susceptibility. 

We can now straightforwardly apply the general for- 
mula (|79p . In spatial Fourier space, this equation applies 
for each wave- vector k. The result is 



Q"(a;,fc) 



1 



M(fe) 



{n{0,k) +iun{Lo,k)) (94) 



Xs{k) + iuj— — — — 
—lu) + Dk 

{Dk^Y + iuDk^ 



{Dk^) 



2^2 



Xs{k). 



(95) 
(96) 



In particular, the imaginary part yields the spectral func- 
tion. 



p""(c.,fc) _ xs{k) 



Dk^ 



TT u;2 + [Dk^y 



(97) 



This equation describes how the transport properties are 
encoded in the low-frequency and low-momentum part of 
the spectral function. The Kubo formula is usually ex- 
pressed in terms of the longitudinal part of the current 
correlator. The latter is related to the density correlator 
by the conservation equation, pl('^, k) = ^p^"{uj, k) (see 
Eq.Hni), 

PL{uJ,k) _ Xs{k) Duj'^ 



UJ TT uj^ + {Dk^y'' 

implying in particular the Kubo formula 



Dy'^ — TT lim lim 



PLiuJ,k) 



(98) 



(99) 



Here ^ 1^ I <ixn{t,x)n{0) is the particle number sus- 
ceptibility. 



3.2.1 Electric conductivity 

The correlators of the electromagnetic current are of par- 
ticular phenomenological interest. We will therefore study 
the operator 



(100) 



In QCD, the explicit expression of the electromagnetic 
current is 



(101) 



In the long-wavelength limit the current-current correla- 
tor is given by the diffusion equation, as described in the 
previous section, however here we wish to include the ef- 
fects of coupling the current to an external electromag- 
netic field, {E,B), and to promote the treatment to sec- 
ond order accuracy - we follow the treatment of [S?' . For 



a linearized theory invariant under parity, Fick's law gen- 
eralizes to 



j = ~DVn + (jE- {(TTj)dtE + kbV x B. 



(102) 



However, a perturbation of the form p,{x)+Ao(x) = does 
not affect the system at all, since both fj, and couple to 
the charge density in the Hamiltonian. Therefore only the 
gradient of this combination can appear in the constitutive 
equation (|102p . and we conclude 



XsD 



(103) 



The conservation law dtJi + V ■ j ~ can be solved for 
n(a;, fe) in the presence of a sinusoidal electric field. Teaney 
and Hong 54 find (for k = (0,0, k)) 



G^^^{io,k)^Dxs 



-ioj + Dk'^ 



Cj^iuj, k) = iDxsUJ - DxsTjuj'^ + KBk^ 



(104) 
(105) 



From here, one can write out Kubo formulae for the second- 
order coefficients tj and kb, 

DxsTj = -^Q^(w, k)^=k=o , (106) 
d 

= gj^G^j^iu, k)^=k=t) ■ (107) 

The coefficient kb can thus be extracted directly by dif- 
ferentiating the static Euclidean correlator with respect to 
fc at cj = 0. An interesting fact about is that it van- 
ishes in a theory whose transport properties are described 
by the Boltzmann equation [511. Thus measuring a non- 
vanishing K B would exhibit a shortcoming of kinetic the- 
ory. The parameter is dimensionless and based on the 
constitutive equation for the current, it describes how a 
flux of particles / der j appears through a surface around 
which there is a non-vanishing circulation of a magnetic 
fleld, / d£ ■ B. Finally, using the conservation law, the 
spectral function for the charge density n is modified from 
Eq. dSIl) to 



p""(w,fc) 



Xs Dk^{l + TjDk^ 



{Dk 



2^2 



(108) 



We conclude this section by noting that the vector 
spectral function determines the emission rate of photons 
by the thermal medium. Denoting the spectral function 
of components /i and v of the electromagnetic current by 
p^v, the real photon emissivity of the quark-gluon plasma 
at temperature T is related to the vector spectral function 
through (see [5S] and Refs. therein) 



P^ui^,k,T) 



(109) 



--\k\ 



where Q{ are the fractional charges of the quarks. Thus 
the conductivity a determines in particular the emission 
of soft photons. Similarly, the thermal production rate of 
dilepton pairs of invariant mass M"^ = a;2 — fe^ reads 



dA^, 



dwdfc-" 



LO^-k^ 



(110) 
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3.2.2 Kubo formulae for the transport of energy and 
momentum 

In hydrodynamics the transport of energy and momentum 
is determined by the correlation functions of the energy- 



momentum tensor T^i,. It is a symmetric, T^^, 



T^^,, and 



conserved operator, 9^7),,^ = 0. The constituent equation, 
which is the analogue of Pick's law (|83)). reads 



T'"' (e + P)m^u'^ + Pg'"' 



(111) 



where u^u^^ 



— 1. The energy- momentum conservation 
equations, together with Eq. 



Ijllip . determine the evo- 
lution of the velocity and temperature fields, u^{x) and 
T{x). To study the linearized modes of these equations, 
one writes it_L, u\\ and e in Fourier components as in (|85|) . 
and decomposes the velocity field into a longitudinal and 
a transverse part with respect to fc, M|| and u±. The re- 
tarded correlator for the momentum density, tt = (e+p)u 
is then obtained along the same lines as Eq. ([M)) . We re- 
fer to appendix C of ^6] for the details and quote the 
leading-order result 



e+p 



j2 - (c,fc)2 + iLuPsk'' 



(112) 



(113) 



respectively for the transverse and longitudinal momen- 
tum density, where 



e + p 



(114) 



determines the damping of sound waves. The spectral func- 
tional is obtained by taking the imaginary part, 

(115) 



1 



rjk 



^ \e+p' 



iuj,k) 



P 



uj'^r.k^ 



(0.2 



(C^nfc2)2 



(116) 



The first line corresponds to the shear channel, the second 
to the sound channel. The former admits a diffusive pole, 
similar to the pole appearing in the vector current correla- 
tor. The latter contains a pole at a; = ztCgk — ^iPk^, which 
corresponds to the propagation of sound waves with ve- 
locity Cs with an attenuation cx e~2^s'= * for a plane wave 
of wavelength A = The sound channel spectral func- 
tion is displayed for two different values of the wavelength 
in Fig. ([5]). Using the momentum conservation equations 
(see [73l) and denoting the spectral function of T^j, and 
Tpa by pf^'^'P'^ ^ one then obtains the Kubo formulae, for 
k = (0,0, fc). 



r](T) = TT lim 



pi3,i3(^^0,T) 





an ^ : t 



lim 



OJ 

p"'JJ(a;,0,T) 
ui 



(117) 
(118) 



Hydrodynamics spectral function Pgo oo^'*') 



a. 

0) 5 



3 4 

'a. 

^ 3 
2 
1 





q / (2;iT) = 0.3 


• 




A 














.1.5 









0.2 0.4 



0.6 0.8 

CO / ( 27tT ; 



1.2 1.4 



Fig. 8. The sound channel spectral function: correlator of Too 
at finite spatial momentum q. This is the functional form pre- 
dicted by hydrodynamics at small momenta and frequencies, 
for Cs = i and TPs = 3^. For comparison with Fig. H15|) . 
the dashed curves are displayed, although they correspond to 
momenta well beyond the validity of hydrodynamics. 

In numerical practice, it is convenient if p{uj)/uj behaves 
smoothly at w — >■ 0. It is therefore important to study the 
spectral function in detail in this region [571155] . One finds 
that 



1 p.ujjiuj,0,T) C 



9 ^ 



- + {e + p)ci5{u) + . 

TT 



(119) 



contains a delta function at the origin (see section [577]) . 
When studying the bulk (scalar) channel it is then prefer- 
able to work instead with the operator 



T 



kk 



(120) 



written here in Euclidean metric, whose spectral function 
at fe = is identical to the spectral function of T*'*' , except 
for not having the delta function at the origin. 

The most general constitutive equation for T'^" has 
been worked out to second order in gradients of the ve- 
locity fields for a conformal theory [53], and more re- 
cently even in the non-conformal case 1471 . In the confor- 
mal case, there are five new transport coefhcients, denoted 
T77, K, Ai, A2, A3. These second-order transport coefficients 
can be obtained from Kubo formulae. 



■qTn 



92 



12.12 



2 9a;2^« 



Gi^'"^^(w, k)uj=k=a 
{cj,k) 



=fc=0- 



(121) 
(122) 



still with the convention k = (0, 0, fc). A physically impor- 
tant effect is that at second order, the speed of sound re- 
ceives a wavelength-dependent correction [53], w = Cs{q)q 
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with 



Csiq) 



1 



rn 



H 

4c2 



(123) 



This provides a way to determine the coefficient Tn if 
the leading-order speed of sound cj. — ^ and Fg are 
known. The coefficient k can be extracted directly from 
the uj = Euclidean correlatofl The coefficients Ai.2,3 
influence the hydrodynamic evolution only at quadratic 
order in the amplitude of the perturbation. Kubo formu- 
lae have been derived for these second-order coefficients as 
well |60j . They involve three-point functions of the energy- 
momentum tensor. Among them we want to point out the 
expression for A3, which like k can be extracted directly 
from a Euclidean correlator, 



A3 = 6 lim 



p,g-i-o dqdk 



^11.01.01/ ,\ 



(124) 



where q and k are aligned along the 3-direction. Both k 
and A3 have a finite limit when the coupling constant goes 
to zero 59; 60]. 

The physical interpretation of the second order coeffi- 
cients is not always straightforward. The coefficients k and 
A3 have been interpreted as a thermodynamic response to 
vorticity in the fluid velocity field [61,60^. The coefficient 
T]j has the unit of time and from its contribution to the 
constituent equation of T^j/, it implies a retardation in 
the response of the fluid to an change in the background 
metric. A relaxation time tji in the sense of the Israel- 
Stewart formalism is reflected in a pole in the retarded 
correlator at Wpoic ~ ^^'''r^- Whether the retarded cor- 
relator actually admits such a pole cannot be answered 
within hydrodynamics, instead it must be answered at a 
more microscopic level, see the discussion in [05] and Refs. 
therein. 



3.3 Quasiparticles and the Boltzmann Equation 

In addition to the pole structure of certain spectral func- 
tions predicted by hydrodynamics, additional character- 
istic features show up if the conserved charges are trans- 
ported by quasiparticles whose mean free path is long com- 
pared to the thermal scale. This is the subject of this sec- 
tion, where we first consider the technically simpler case of 
a diffusing heavy quark, and then discuss the transport of 
quark number, energy and momentum by the constituents 
of the quark-gluon plasma. 



3.3.1 Diffusion of a heavy quark 

The diffusion of a heavy quark (one has in mind the charm 
or bottom) in the quark-gluon plasma is characterized by 



* See Eq. (|49p . where in view of p88|l no subtlety occurs 
when taking the limit oj — >■ 0. Note however that a quartic 
divergence must be subtracted from G^^'^^. 



a time scale M/T"^ which is long compared to the ther- 
mal time scale of 1/T. For this reason it is widely be- 
lieved that a classical Langevin equation should appro- 
priately describe the thermalization of heavy quarks [35] . 
The equations of motion for the latter are 



dx p 

H ^ M 



^^m-wit), (125) 
t'). (126) 



For a given £_{t), the equation is easily solved to give 



Then 



p(0=e-''*[p(0)-K/ods|(s)e''^ 



lim {p,{t)pj{t)) = —S,j. 

t— Voo Z71 



(127) 



(128) 



The equipartition of energy requires to be |T in equi- 
librium. The drag and fluctuation coefficients are then re- 
lated by the fluctuation-dissipation relation (established 
by Einstein in 1905), 



K 



V ■ 



2MT 



(129) 



Furthermore, the mean square distance is also easily worked 
out. 



Ux^it)) 



t- 



i(l-e- 



■Tit 



)-i(l 



-r;^^2 



(1- e 



(130) 



For large times, the first term dominates. Recall that the 

diffusion equation ([57]) yields a density of particles n{t, k) = 

■2 2 
n(0, fc)e~^'° *, or n{t,x) cx e~TOt for an initial distribu- 
tion localized at the origin. Therefore the mean square 
radius of the particles after time t is ■^{x'^) — 2Dt, and 
comparing with p30p , one finds that the Langevin process 
leads at late times to a diffusion with diffusion coefficient 



D 



T 

Mr]' 



(131) 



For a thermal initial distribution of momenta, {pi{0)pj (0)) 
MTSij, expression (|130p simplifies to 



x\t))=2D[t-lil-e-^')] 



(132) 



This equation describes both the early-time directed mo- 
tion, i(a;2(t)) = -jw^i^, ^v"^ = -p., and the late-time diffu- 
sive motion, l{x^{t)) = 2Dt pi. 

We follow the treatment of |63j to establish the con- 
nection between the distribution of heavy quarks after a 
time t and the retarded correlator of the heavy quark cur- 
rent. Let P{t, x) be the probability that a heavy quark 
starts at the origin at t = and moves a distance x over 
a time t. If the distribution of heavy quarks at time zero 
is iV(0, a;), at time t it will be given by the convolution 



iV,.,x)^/dx'P(,,x-x')«(0,.0, 



(133) 
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or equivalently 

N{t,k) = P{t,k)N{0,k). 
Applying the general rule (l79l) . we thus find 



1 + iuj 



(134) 



(135) 



If one assumes the noise to be Gaussian distributed, then 
the probability distribution P(i, x) is Gaussian 63 , with 
a width given by Eq. (jl32p . and therefore so is P{t, k). The 
static susceptibility Xs{k) can be obtained by noting that 
the initial phase-space distribution of the heavy quarks 
is described by f{x,p,t = 0) = e'3(A'(a')-A^-p'/2M) rp^_ 
king into account both quarks and antiquarks, the static 
susceptibility is independent of k to linear order in the 
perturbation ^(x) = /io + S^i{x), 



Xs(fe) = - 



MT 



2tt 



3/2 



e"^*^cosh;3/zo, 



(136) 



where u = ANc is a multiplicity factor. Eq. p35p and <\VdQ\ 
thus provide all the ingredients to compute the kinetic 
theory prediction of Gn at low frequencies and momenta. 
Fig. from [53], displays the spectral function of the 

longitudinal current correlator {pL{uj,k) = ^^ImG^"(a;,fc 
denoted by pjj on the figure). In particular, the k = 
spectral function takes the form 



V 



(137) 
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Fig. 9. Heavy-quark spectral function of the longitudinal cur- 
rent correlator as predicted by the Langevin equation [63] as 
a function of a scaled frequency Co = ujD{M/T) for various 
values of a scaled momentum k = kDM/T. The solid lines are 
the predictions of the Langevin equations, while the dotted 
lines are the predictions of the free theory. The dash-dotted 
line shows the fc = result of the Langevin equation. 



By comparison with the hydrodynamic prediction (|98p . 
we see that the kinetic treatment predicts not just the 
intercept of pl/'^i but also its local functional form, a 
Lorentzian. The width has the interpretation of an inverse 
relaxation time. It is a long-known fact that kinetic theory 
establishes a connection between the diffusion coefficient 
D — -p^ and the relaxation time 77"^. Already Drude's 

classical description of the electric conductivity in metals 

2 

around 1900 showed that a/{Tn) = where n is the 
density of electrons and m their mass. 

The spectral structure that is obtained from the Langevin 
equation is expected to arise for a sufficiently heavy dif- 
fusing particle. Conversely, the presence of a transport 
peak allows one to define the quantities appearing in the 
Langevin equation directly from the spectral function [64 . 
The effective mean-square velocity is then given by 



If one inserts the spectral function (jl37p derived from the 
Langevin equation, one recovers Einstein's relation (|129p . 
A weak-coupling calculation shows that while k{M) and 
D are only weakly dependent on Af, the drag coefficient 
•q ~ T"^ /Mx a power of the coupling constant is paramet- 
rically small compared to the medium time-scale. This jus- 
tifies a posteriori the assumption that the transport peak 
is narrow. The momentum diffusion coefficient k, converges 
to a finite value when the static limit M — > 00 is taken. 
Note that although the peak in pl{^) becomes arbitrarily 
narrow in that limit, the limit M — >■ 00 cannot be inter- 
changed with the limit w in the definition of n. 



= — / Pi w , 

Xs J^A ^ 



(138) 



where A is a cutoff that separates the scale 77 from the 
correlation time of the medium (which is typically of or- 
der T, or gT at weak coupling). The 'kinetic mass' Mkin 
is further defined so as to satisfy the equipartition the- 
orem, Mkin('U^) = 3T. Finally, the momentum diffusion 
coefficient k{M) can be defined as 



The formula (|139p can be understood heuristically as 
follows : PL is the spectral function for the operator jk 
at vanishing spatial momentum. Classically, if the quark 
is heavy, this operator measures a single quark's velocity, 
Vk — J dx jk ■ Therefore the operator M J dx mea- 
sures the force acting on the quark. We now recall that k 
measures the size of the force exerted by the medium on 
the quark, Eq. (|126p . We then have to evaluate the clas- 
sical force-force correlator, whose quantum analog is the 
symmetrized correlator. Combining these ingredients, one 
reaches 



K — — lim 

6 ^LU^O 



k=l 



]\/r2 poo 

lim I dtt 



L)(t-t') 



(140) 



k{M) 



— wpl(w) 



i)<|w|<yi 



(139) 



/ r djk{t,x) djfc(i',0),> 



dt' 
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A second, rigorous derivation was given for this formula 
in [64]. Through the fluctuation-dissipation formula (p9|) . 
the symmetrized correlator can be expressed in terms of 
the corresponding spectral function, which then leads back 
to formula (jl39p in the heavy-quark limit. As one might 
expect based on the classical argument laid out above, the 
correlator (jl40l) can be expressed in the Heavy-Quark Ef- 
fective Theory (HQET) as a force-force correlator, where 
the leading force is the chromo-electric force gE, 

K=-V lim — / dtdx(U(f)''gE''$-§''gE''e)(t,x), 

{4>^gE'^$-0^gE'^e){O,O)}). (141) 

where 9 and are two-component spinors of HQET and 
the w — >■ limit has now been taken. The corresponding 
Euclidean correlator reads, after evaluating the fermion 
line contractions, 



^^^^ _ ^ReTT{Uil3,t)gEkit,0)Uit,0)gEkiO,0))^ 

~ -3 (ReTr[/(/3,0)) ' ' 

(142) 

where the color parallel transporters U{t2,ti) in the fun- 
damental representation are propagators of static quarks. 
In particular the denominator of (I142p is the Polyakov 
loop. The momentum diffusion coefficient is given by the 
low-frequency limit of the corresponding spectral function 
via Eq. 



,. 2ttT 
lim / 



,HQET 



(143) 



The obvious advantage of this formulation is that the large 
scale M has disappeared from the problem. The spec- 
tral function pHQET j-^g^g ]-,ggj^ studied in detail at next- 
to-leading order in perturbation theory [SS]. Remarkably, 
even in the weak-coupling limit, the function is smooth 
as small frequencies. This is in contrast with the narrow 
transport peaks that are found at weak-coupling in e.g. the 
shear channel. This property is a clear advantage for nu- 
merical studies of the spectral function, as we will see in 
section 



3.3.2 Kinetic theory in other channels 

When studying the diffusion of a light quark number, or 
the transport of transverse momentum, the applicability 
of kinetic theory is not a priori guaranteed. It is only valid 
if the medium admits quasiparticles that are sufficiently 
long-lived. In the heavy-quark diffusion case, the existence 
of a quasiparticle is guaranteed by the separation of scales 
M/T. Typically, although perhaps not necessarily, quasi- 
particles arise in weakly coupled systems. In 0^ theory, the 
quasiparticle carries the same quantum numbers as the 
particles at zero-temperatures. At sufficiently high tem- 
peratures in QCD, due to asymptotic freedom, quasipar- 
ticles are expected to emerge, with the quantum numbers 
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Fig. 10. Kinetic theory prediction for the vector current cor- 
relator in quenched QCD for massless quarks f54]. The mo- 
mentum is directed in the z direction, fc = (0, 0, fc), and 
fiF ~ g^Cpmjj log(r/m_D)/87r is the drag coefficient of a quark 
in leading-log approximation. The thick dashed curve corre- 
sponds to k = 0. The thin dashed lines show the result of the 
free Boltzmann equation [661154) . Note that the normalization 
of p differs from our convention by a factor 2n. 



of quarks and gluons. At low temperatures, well below the 
crossover region, we expect the well-known hadron reso- 
nances to be good quasiparticles. In that regime the pa- 
rameter 1 /Nc can be used to provide an expansion around 
the non-interacting hadron gas {Nc is the number of col- 
ors). 

The physical picture one has in mind is that of a par- 
ticle propagating freely until it makes a 'collision' with 
a constituent of the medium, upon which its momentum 
or even its nature may be changed. The collision gener- 
ally has to be described by a quantum-mechanical matrix 
element. The propagation of the quasiparticle between col- 
lisions can be treated as classical, provided the time until 
the next collision is long enough that any coherence effects 
can be neglected. The evolution of the quasiparticle can 
then be written as a convolution of transition probabil- 
ities, rather than amplitudes. In this approximation the 
evolution is described by the Boltzmann equation for the 
probability distribution f{x,p,t). 
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To ilUustrate these general remarks, we review a recent 
kinetic theory treatment of the diffusion of light quarks |54) . 
To treat this problem one needs the Boltzmann equation 
in the presence of an external field. Recalling that the 
Lorentz force on a charged particle is = QF^^v^^, the 
Boltzmann equation (for a single flavor charged under the 
externally applied gauge field) reads 



QF^''p,-^^r = C'^[f,p]- (144) 

Here Q denotes the electric charge of the quark and below 
we will employ the notation 



n^(w) : 
n^(w) 



1 



1 



(145) 
(146) 



often with the understanding of an ultrarelativistic dis- 
persion relation, i.e. = n^{uj = |p|). Specializing to a 
spatial gauge potential, = (0, A), and linearizing the 
Boltzmann equation yields in Fourier space (cj, k) 



{-iu + iVp ■ k)5f°-{uj, k) 



(147) 



-i^n{l-n)QA,^^C''[Sf,p\. 

hirt 



Since the coupling of the gauge field is linear in the charge, 
it only influences the difference of quarks and antiquarks. 



-lUJ 



+ iVp ■ k)5f 
-iujn{l — n)2QA 



{LO,k) 



(148) 



At this point one needs to specify the collision kernel C^'' . 
The quarks and gluons contribute a 'loss' term (1st line) 
and a 'gain' term P^I54j . 
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_2^<(l + n.^ 
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dfc 



P J (27r)3 



5/(p)=n^(l-n^)x(p), 
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(150) 
(151) 

(152) 



Hong and Teaney solved the Boltzmann equation (I148P 
by expanding functions of momenta in real-valued linear 
combinations of spherical harmonics. The radial part of 
the momenta is discretized, and the problem is reduced to 
inverting a large matrix. The out-of-equilibrium expecta- 
tion value of the current is determined from 5f via 



{Ji)A = QVs 



dp Pi 
(27r)3 Ep 



5f 



(153) 
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Fig. 11. Spectral function of T^y in the shear channel, as pre- 
dicted by the Boltzmann equation [54] in leading log approx- 
imation. The thick dashed curve corresponds to fc = 0, and 
= (;^Ca?w|) log(T/m_D)/87r. The thin dashed curves rep- 
resent the free theory result |67ll54j . The normalization of p 
differs from our convention by a factor 2-k. 

and from there the retarded correlator is obtained via the 
frequency version of the linear response formula (|75p . 



(J^(^,fc))A 



G^"(tj,fc)A,(w,fc). 



(154) 



The result is shown in Fig. ((10)) . In the rescaled variables, 
the functional form of the curves corresponding to dif- 
ferent momenta are similar to those obtained for heavy 
quarks. Fig. 

A similar approach was used in [51] to determine the 
correlators of T^^. Now the perturbing field is gravita- 
tional rather than electromagnetic, and consequently the 
Boltzmann equation has to be generalized to non-fiat back- 
ground metrics. Another important aspect of the calcula- 
tion is the selection of appropriate boundary conditions. 
Hong and Teaney argue that the appropriate boundary 
condition for x(p) is a Dirichlet boundary condition at 
p = 0. The reason is that the rate of soft gluon emis- 
sion is parametrically high compared to the transport time 
scales. At high momenta the solution that grows exponen- 
tially is discarded by implementing the appropriate finite- 
difference scheme at the highest discretized momentum. 
The result for the shear channel is shown in Fig. (|TT|). By 
the Kubo formula (jll7p . the shear viscosity in leading log- 
arithmic approximation can be read off from the intercept 
of the fc = curve (parametrically, rj/ s ^ [g'^ \og{l / g)]~^) . 



3.4 Leading order results for the shear and bulk 
viscosities 

The conductivity, as well as the shear and bulk viscosi- 
ties, have been obtained beyond the leading-logarithmic 
order, namely to full leading order in the gauge coupling 
g [T51I68] . The method follows the steps sketched in the 
previous section, namely a Boltzmann equation is writ- 
ten down and subsequently linearized. In addition to the 
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ingredients necessary for a leading-logarithmic accuracy, 
the collision kernel must be treated more carefully at low 
momenta, and the Landau-Pomeranchuk-Migdal (LPM) 
effect, which involves collinear emissions and their destruc- 
tive interferences, must be taken into account |51) . 

Arnold, Moore and Yaffe (AMY) solved these linear 
integral equations using a variational approach (a method 
introduced earlier in ,69)). Numerical results for the shear 
and bulk viscosity are displayed in Fig. The results 

are divided by the entropy density, effectively normalizing 
them by the number of degrees of freedom. The ratio 77/s is 
about 30% lower in the pure gauge theory than in iVf = 3 
QCD, primarily because the larger color charge carried by 
the gluons makes them equilibrate faster. 

Numerically, r]/s is about 2.0 for — 0.15, but the 
accuracy of the calculation deteriorates for larger values 
of the coupling. In the weak coupling regime, the bulk 
viscosity, of order 5^ in the gauge coupling, turns out to be 
negligible compared to the shear viscosity, see the bottom 
panel of Fig. ([H]) . 

Recently, the treatment of the inelastic process gg ^ 
ggg in the shear viscosity calculation has been scrutinized 
by Xu and Greiner (XG) [7(J|. They showed that at mod- 
erate values of ag this process can numerically dominate 
the elastic gg O gg scattering due to soft, not neces- 
sarily collinear emissions. The calculation was revisited 
very recently by Chen et al. j71] without relying on the 
soft gluon bremsstrahlung approximation. The inclusion 
of gluon emissions appears to lower somewhat the ratio 
77/s, but a full next-to-leading order calculation has not 
been performed yet in a non-Abelian gauge theory. 



3.5 Perturbative predictions for spectral functions and 
operator production expansion 

We begin with the energy-momentum tensor correlators. 
At frequencies a; 3> T, straightforward perturbation the- 
ory is expected to apply. Here we give some results for 
the SU(A^c) gauge theory in the limit of zero coupling. 
In the non-interacting limit, gluodynamics is of course a 
scale-invariant theory. There are thus three independent 
correlators of the stress-energy tensor, which were calcu- 
lated in ^7] . To describe them, we assume in the following 
expressions that w, q > 0, g = (0, 0, q). We define the func- 
tional 

, Piz) sinh ^ , , 



dz 



2T ^w-^ii 2T 

P{z) sinh §p 
cosh ||t — cosh 2y 



Thus 1\P\ is a function of {uj,q,T). Then the spectral 
functions read, respectively in the shear channel, the sound 
channel and the tensor channel. 
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Fig. 12. The shear (77) and bulk {(,) viscosities scaled by the 
entropy density, in full leading order of perturbation theory 
([5T][T5] and [BS] respectively). 



where (Ia = — 1. In addition, the correlator of the trace 
anomaly reads 



pe,e[uj,q,l) -[—^) Tu^y"^ "9) 



(159) 



^ 2T , sinh(w + g)/4r 

-e{q -uj) + — log — ^ ^-f— 

q q smh \u) — q\/4:l 



For P a polynomial, the quantity I[P] can be expressed in 
terms of polylogarithms i67j. The correlators of T^i, can 
be expressed in terms of polylogarithms up to order 5. The 
spectral functions are displayed for q = ttT in Fig. (fT3|) . 
Using expressions of this kind, the second-order transport 
coefficient k can be extracted via the Kubo formula (jl2ip , 
yielding [5 7) , after subtraction of the vacuum contribution. 



^^2 ^ 5s 
18 87r2T' 



(160) 



Even in the non-interacting limit, the functional form 
of the spectral functions is relatively complicated, due to 
the three scales uj,q,T involved in the problem. The ex- 
pressions simplify significantly when the spatial momen- 
tum q is set to zero, for instance (without summation over 
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Free spectral functions p^vpo in SU(N) gauge theory 




x[6ji/(11a3 N)]'^ nn,vv 
x(1/4) 00,00 
12,12 
13,13 



q=(0,0,7iT) 



10 15 
co/T 



20 



25 



Fig. 13. Non-interacting spectral functions in the SU(A'') 
gauge theory. More precisely, q = ttT and the functions plotted 
are ^^rpi \p{'-^i 1i T) / tanh(a;/2r) — p{u), q, 0)] for four different 
channels. The sound channel poo,oo admits a discontinuity at 
threshold, uj — q. 



the indices fi, v) 



Any 



A 



12.12 — ^13,13 



tanh ^ujP 

1 dA 
10 (47r)2' 

^12, 12 ~ -Bi3,i3 = j ^a; 



— + B^,r^LjS{Lu), (161) 



Ag = dA 

Be = 0. 



lla^iV 
3(47r)2 



(162) 

(163) 

(164) 
(165) 



The 5-function at the origin corresponds to the fact that 
gluons are asymptotic states in the free theory and implies 
an infinite viscosity. 

Next-to-leading order corrections are known in some 
cases. In the bulk channel p^^^^u, they have been cal- 
culated recently [72j, see Eq. (jl8ip . In the shear chan- 
nel, the term in the large-frequency behavior of the 
spectral function ^13^13(0;, 0, T) is corrected by the factor 

[1-^] mm ' 



Expressions similar to Eq. (jl56p were obtained earlier 
by Aarts and Martinez Resco for the free spectral func- 
tions of dimension-three quark bilinears [6B] in terms of 
polylogarithms. In this case polylogarithms up to order 3 
appear. At g = 0, the expressions again simplify consid- 
erably. For the vector current, the result is for one flavor 



(see for instance [751l63j ) 

3 



(166) 
(167) 



7^^(^-2™)\/l-( — ) 



47r2 

+ 



2m^2\ 9 w 
UJ" tanh — -. 
AT 



The susceptibility and average velocities have simple ex- 
pression in the massless and in the heavy-quark limits. 



J.2 



4Na 

T 



2ir 



1 

-pm 3T 



TO = 



(168) 



At high frequencies lo ^ T, the leading power in uj 
of the spectral function gets modified by radiative cor- 
rections. In this frequency regime, they are independent 
of the temperature and can be calculated in the vacuum. 
Famously, the vector-current spectral function pa is then 
modified by a factor (1 -I- ^) [32] ■ The quark number sus- 
ceptibility Xs also receives calculable corrections at high 
temperature, Xs = ^^^(l — ^^^^) in the massless case 
[75] . The broadening of the delta function in pa into a 
peak of width on the order of the inverse relaxation time 
however requires either an explicit resummation procedure 
or a kinetic theory treatment, see section p. 3D . 



3.6 Thermal Correlators from the Gauge/Gravity 
Correspondence 

In this section we summarize some of the results on spec- 
tral functions that have been obtained in the strong cou- 
pling regime using the holographic principle. We will re- 
strict ourselves to presenting some of the spectral func- 
tions of the A/" = 4 SYM theory. Our main goal here is to 
contrast their analytic structure with the analytic struc- 
ture found at weak coupling using kinetic theory. For the 
interested reader, the method of computing spectral func- 
tions based on the holographic principle [771l78[l79l[5ni[5S] 
has been reviewed in [8ni49l[50| . 

The spectral function of the R-charge correlator at 
vanishing spatial momentum, remarkably, can be written 
in terms of elementary functions |82j . 



sinh 2y 



27rry cosh 



2T 



cos 



(169) 



2T 



The static susceptibility is given by Xs — 
diffusion constant D = 

27r_/ 



g , and the 
can be read off from this for- 



mula by using the Kubo formula, Eq. This function 
has a series of poles along the diagonals in the lower half 
of the oj complex plane. 



Wo 



n ■ 27rr(Tl 



(170) 
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Fig. 14. Spectral function of shear stress at strong coupling in 
the A/" = 4 SYM theory at vanishing spatial momentum. The 
dashed curve shows the zero temperature result (|171|l normal- 
ized by the same factors. Due to a non-renormalization theo- 
rem in these channels, the zero temperature spectral densities 
in the free and interacting theories are equal |83II84) . At fi- 
nite temperature the kinetic theory peak observed in Fig. Hill 
fc = curve) does not exist in the strongly interacting theory. 
Reproduced from [56] . 



The shear channel spectral function {p{ui)/uj) displayed 
in Fig. [56,80 smoothly interpolates between the vac- 
uum spectral function 



P?2l2(^,'7 = 0,T = 0) 



647r2 



(171) 



at large frequencies and a finite value at w = 0, which 
yields the shear viscosity rj = ^N^T^ via the Kubo for- 
mula (I117p . Combining this with the entropy density of 

s = ^^N^T^, one obtains the ratio rj/s ^ l/in (see Eq. ([3]) 
in the Introduction). 

Numerically no particular structure or excitation is 
visible in the spectral function at any frequency. This is 
in stark contrast with the peak structure predicted at the 
origin by the Boltzmann equation, see Fig. ([TT|). One con- 
cludes from figure (|T4)) that the conformal plasma of the 
J\f = A SYM theory does not admit quasiparticles. There 
is no separation between the transport time scale tjj and 
the thermal scale ttT. In fact, the second-order transport 
coefficients were calculated in the strongly coupled N = A 
SYM theory ([59], see also [85] for r/j), with the result 



Tn 



log 2 



2ttT 
ttT 



(172) 
(173) 



The sound channel spectral function, namely the spec- 
tral function for the energy density Tqq at finite spatial 
momentum, is diplayed in Fig. (|15p. Comparison with 
Fig. (jSj) shows that in this strongly coupled theory, hydro- 




Fig. 15. Spectral function of the energy density at spatial 
momenta obtained from the AdS / CFT correspondence for the 
— 4 SYM theory. The quantity displayed is Tff^ in infinite- 
coupling, large- A''c SYM theory, computed by AdS / CFT meth- 
ods. The X-axis variable is w= u}/2tiT and the curves cor- 
respond to q/2-KT = 0.3,0.6,1.0 and 1.5. Figure reproduced 
from EOl. 



dynamics provides a good description of the spectral func- 
tion up to momenta as large as ttT. In a weakly-coupled 
theory, this range of validity is parametrically smaller. 

Assuming the holographic principle to hold for other 
theories as well, non-conformal theories can also be stud- 
ied by these methods. In the past few years, a program has 
been undertaken to tune the properties of the 'bulk' theory 
to reproduce the thermodynamic properties of QCD on 
the 'boundary', and then predict those properties such as 
transport coefficients that are difficult to extract from first 
principles [S51I871I88] . These models can give qualitative 
insight into whether the shape of the spectral functions 
displayed in Fig. ([Ti)) and are robust again modest 
modifications of the theory. A longer-term goal could be 
to understand the crossover between the functional form 
observed at weak coupling and that found at strong cou- 
pling. 

We remark that in the 'AdS/QCD' models the shear 
viscosity to entropy density ratio remains 1 /47r at all tem- 
peratures, which shows that the breaking of conformal in- 
variance is not of the most general form. A bound that 
is found to be observed by these theories [SS] is Buchel's 
conjectured bound [55] . 



> 2 



(174) 



It implies that the bulk viscosity must become at least 
comparable to the shear viscosity when the speed of sound 
becomes small, as happens for instance at certain phase 
transitions. 



3.7 Thermal Sum Rules 

We briefly remind the reader of the derivation of sum rules. 
The first step is to derive a dispersion relation [SD]. For 
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this purpose it is useful to introduce the function 



Gfl(w) EE die*"* G{t) Ini(cj) > 
Ga{oj) = - dt e»"* G{t) Im(w) < 

(175) 

It is analytic in w everywhere except along the real axis. 
We now apply Cauchy's theorem for a contour C which is 
the sum of two large semi-circles, as depicted on Fig. ([T5|) . 
In order for the contributions at infinity to vanish, the the- 
orem is applied to Q{uj) / P{ijj), where P{uj) is a polymonial 
of sufficiently high order n with simple zeros at a; = uj^. 
The contour integral then yields 

1 / G{z) dz _ g(a;) x - G{i-Ov) 

{z-uj)P{z) ~ P{uj) ^ {uj^-uj)P'{uj^y 

(176) 

By construction, the contribution of the semi-circles van- 
ishes in the limit of infinite radius. Now recalling the def- 
inition of the spectral function, the left-hand side is equal 
to 

dz p{z) 



{z-Lo)p{zy 



(177) 



Therefore, noticing that t/(w) = Gii{uj) for Im(uj) > 0, we 
obtain 

Gfi(a.) = Pico) r , "^^P^l ^ + g(a.), (178) 



where 



-P(u) 



{uji, - uj)P'{lu^) ' 



(179) 



is a polynomial of degree n — 1. In practice it means that 
n constants have to be calculated before formula (|178l) 
determines the retarded correlator For instance, for the 
vector correlator, Eq. (|lip . one subtraction is sufficient 
to make the dispersion integral converge. We reserve the 
name sum rule to the case where Gii{uj) can be calculated 
in an independent way, for instance in terms of thermo- 
dynamic potentials. Eq. (jl78l) can then be interpreted as 
a global constraint on the spectral function p(w), which 
may help to parametrize it in a consistent way. 

Since the derivation of a sum rule involves integrating 
the spectral function over all frequencies, it is necessary to 
know its large frequency behavior in order to ensure the 
integral's UV convergence. This leads us to the subject of 
the operator product expansion (OPE). In thermal field 
theory, the spectral functions typically rise at large fre- 
quencies with a positive power of frequency. The simplest 
subtraction to implement is then to take the difference 
between two temperatures. This reduces the asymptotic 
power by four units, due to the absence of gauge invari- 
ant operators of lower dimension (chiral symmetry implies 
that the operator ^ijj appears accompanied by a power of 
the quark mass). 

Correlation functions deep in the Euclidean region ad- 
mit an operator-product expansion. Through the disper- 
sion relation, this leads to a large frequency expansion for 
the spectral function [5T]. Leading-order results for the 




Fig. 16. Contour chosen in the derivation of thermal sum 
rules, see Eq. (|178|) . 



OPE of current and energy-momentum tensor correlators 
were obtained in [S2] ■ Recently, a full two-loop calculation 
of the correlator oi 9 = T^^ has been performed [75] in 
Yang-Mills theory (see also the closely related work [55]). 
In addition to giving the radiative corrections to the lead- 
ing power, this calculation provides the Wilson coef- 
ficients of the 9 and 6*00 = Tqq ^ j9 operators in the op- 
erator product expansion of the trace anomaly correlator. 
We quote the uj :s> ttT result. 



4c2 (47r) 



(180) 



2(e-hp)7e;e + p(t^) 
2(e-3p)7e;e-3p(w) + Of ^ 



76(;i(w)=g + . , s n I — In 



(47r)2 y 3 w2 3 
-0(/) , 



7,;, _3p(^) -.9^ + 0(5^), 



4 y I 



6n = 



and bi = 



(181) 

(182) 
(183) 



where ce - - , i^u - -^ii^ <^^^^ - g^^fy 
and is the dimensionless renormalized coupling con- 
stant, evaluated at the MS scheme renormalization scale 
ji {p? = p?e^^ /^tt). Expressions such as Eq. (I18ip are 
very useful to constrain the form of the spectral function 
at high frequencies, particularly since the thermodynamic 
potentials e and p can be evaluated non-perturbatively. In 
particular, if the vacuum correlator is subtracted from the 
thermal correlator, then the leading asymptotic behavior 
is determined by two linear combinations of energy density 
and pressure. This has recently been exploited to constrain 
the spectral function pg [58], at which time however only 
the leading order behavior of the Wilson coefficients was 
known [92] . 
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The calculation shows that (a) the scalar and pseu- 
doscalar correlators do not have contact terms (the cor- 
relators vanish in the limit w — >■ oo) and (b) the integral 
of Ap/uj is UV-finite. Using these facts, one can derive a 
sum rule in the SU(7Vc) Yang-MiUs theory [Mll^ . First, 
the following expression 



.(0) 



d 



3p 



dT 



(184) 



is obtained for the (w, q) = Euclidean correlator, where 
e and p are the energy density and pressure of the finite- 
temperature system. We now convert identity (jl84p into 
a sum rule for the bulk spectral function with the goal to 
constrain the latter. Combining equations (jl84p and (jH), 
we obtain 



— [pe(w,r) - pe(w,0)] 



3p 



rp4 



(185) 



Next, we separate an infrared contribution of the form 
uj6{ui) in pe from the rest of the spectral function [57] . 
Starting from Eq. (jll6p . we recall that for any smooth 
function such that f{x)dx — 1, ^/(x/e) provides a 
representation of the delta function. Here we will exploit 
this fact for the function 



/(^) = 



1 



7r(62 - 1) 



(a;2-62)2+; 



- 1 



(186) 



and b — Cs/iFsk) and e — Fgk^ . In this way one finds that, 
in the sense of distributions. 



lim 

fc->-0 



P33,33('^,fe,7') e + p 



U> TT 

In the shear channel on the other hand, 



r, + {e+p)clS{Lo) + ... (187) 



Pl3,13(w,fc,T) 



V 

TT 

t 

TT 



e+p 



(188) 



21,2 



rj k 

e + p 



and the delta function is suppressed by a power of k . 
Thus combining the shear and sound channels. 



1 Pii,kk{i^,0,T) 
9 



P33,33(t^,0, r) _ 4 pi343(tJ,0,T) 

3 



OJ 



c 



= - + {e+p)ciS{uj) + 



Using the exact relations 
poo,oo(w,0,T) 



Poo,fefc(w,0,T) = -3{e + p)ujS{uj) , 



(189) 

(190) 
(191) 



one finds 



9(0;, o,r) 



OJ 



9C 

TT 



-(3c2-l)25(c.) + 



We introduce the operator 

a = i-^/' / d^'x {T^^{x) + {icl - 1)Tqo{x)) . (192) 



Its spectral function, which we denote by p^,, satisfies 



pe(w,0,r) = p*(w,0,T) 



P 



{icl - if OJ 5{uj) (193) 



and, in view of Eq. (I192p . is free of the delta function at 
the origin. 

We are now ready to turn Eq. p85D into a sum rule for 
p*, the spectral function of the operator Oi, (jl92l) . which 
is smooth at the origin. Using Eq. (|193p . one finds 

— [p,{oj,Typ,{oj,{))] = 3(l-3c2)(e+p)-4(e-3p) . 

OJ 

(194) 

Eq. (jl85p was already obtained in [SS], but the presence 
of the oj5{oj) term in pg was missed. In [571 , a finite spa- 
tial momentum k was used as an infrared regulator and 
the contribution of the zero-frequency delta function cor- 
rectly identified for the first time. If one insists on ex- 
pressing Eq. (I194p in terms of pe, one should strictly write 
lim,_o/_°l:^Z\pe(a.,r) = 3(l-3c2)(e+p)-4(e-3p). 

In the shear channel, a similar sum rule applies in the 
Yang-Mills theory, as first pointed out in [57] . It was later 
also studied in lattice regularization [96] . There it was 
given in the form 

— Z\pi2.i2(w,P = 0,T) = -e(r) (195) 
-00 ^ •J 

- lim AG{oj,T), 



where 



AG{uj,T) = j d^ 



-|Too(x)Too(0) (196) 



+ \{T^i{x) - T22(a:))(rn(0) - T22(0)) 



T-O 



was proved to be finite. The form (|195l) is equivalent to 
the form given in (57| . The sum rule remains incomplete, 
because a finite, but yet undetermined contribution of the 
trace anomaly can affect AG{oj — )■ 00, T). The determina- 
tion of its OPE coefficient requires a two-loop calculation. 
To our knowledge, this contact term can not be deter- 
mined by Ward identities. 

Another important example is the difference of thermal 
vector and axial- vector current correlators, which general- 
izes the Weinberg sum rules of the sixties. The difference 
of the corresponding spectral functions obeys sum rules 
Vfe in the chiral limit [97] . 



Aoj OJ 



^ 



{pl{oj,k)~pi{oj,k))=0. (197) 



diooj {pl{ij,k)~ pi{io,k)) =Q. (198) 
dww (p^(a;,fc) - p^(w,fc)) = 0. (199) 
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A remarkable point is that the difference of the vector 
and axial- vector spectral functions falls off rapidly at large 
frequencies. In the OPE analysis, the correlators first differ 
by a four-quark operator, which means that the difference 
of spectral functions falls off as w"^ at high frequencies. 



3.8 Correlators in Finite-volume 

All the analytic treatments covered in this section are per- 
formed in infinite spatial volume. However, the Monte- 
Carlo calculations reviewed in section ^ are necessarily 
performed in finite volume V = L^. The boundary con- 
ditions are normally taken to be periodic. As is manifest 
from the formal spectral expression (l46l) of the spectral 
density, the latter is a distribution, which consists of a 
collection of delta functions with a weight given by the 
matrix elements of the current under investigation be- 
tween energy eigenstates. The spectral weight at the origin 
oj = is thus determined by transition matrix elements of 
the current between states whose energy differs only in- 
finitesimally. Symmetry considerations thus do not result 
in as stringent restrictions as at zero temperature. For in- 
stance, for any energy eigenstate above the two-particle 
threshold, there typically exist energy eigenstates whose 
energy differs only by an amount of order 1/L and which 
are connected to it by a non- vanishing matrix element. 

Clearly one cannot literally apply the Kubo formula to 
the finite- volume spectral function, since this would given 
a undefined value. The situation is similar in this respect 
to the discussion of the spectrum of the Dirac operator. In 
that context, the spectral density at threshold in infinite 
volume is equal to the chiral condensate —{iptp). In finite 
volume, the Dirac spectrum is discrete and therefore it is 
necessary to send the volume to infinity (before setting 
the quark mass m to zero) in order to extract the value of 
the chiral condensate. However, if the density is given by 



p(A,m) = f]5(A~Afe), 



(200) 



fc=i 



the difficulty can be dealt with by working with the inte- 
grated spectral density |98) . 



i/(M,TO) = y/ dXp{X,m), yl=VM2 



(201) 

A finite spectral density manifests itself by a linearly rising 
i^(M, m) as a function of M on a scale much larger than 
the typical eigenvalue interspacing, but much smaller than 
the scale on which a curvature arises. This behavior is 
readily seen in simulations of QCD in large but finite vol- 
ume [55] • In the case at hand of thermal spectral functions 
of gauge invariant operators, a natural way to deal with 
the singular nature of the spectral functions is to work 
with the coordinate-space retarded correlator G^^{t). For 
B = A = A\ inverting the Fourier transform ([32t yields 



G{t) = 2 r 

JO 



duj sin(cji) p{io). 



(202) 



Consider for illustration the quantum Hamiltonian of 
a free particle coupled to an assembly of harmonic oscil- 
lators {[q,p] = i, [qn^Pm] = iSnm) 99,100;, 



2m ^ -^-^ 



(203) 



The Heisenberg equations of motion Y — i[H,Y] are eas- 
ily obtained (first for the annihilation operators a„ = 

""v'2it^^" ' ^^"^^ '^'-'^ P ^^'^ ^^'^ ^^^^ form 



q{t) = p{t)/m, 

m = m - 



ds fit- s) q(s) , 



2i 



a„(0)e 



4(0)e 



iuJnt 



fit) = ^nlcosujnt. 



(204) 
(205) 
(206) 
(207) 



The operator equations for q and p resemble formally the 
classical Langevin equation with a memory kernel given by 
fit). The commutator-defined correlator of the operator 
^(<), which plays the role of the random force term, is 
given by 

dfit) 



GiHt)^m),m] = - 



dt 



and, based on Eq. (j202p . the spectral function by 

2p(w) 2r/ N 



(208) 



(209) 



It has the expected form for a spectral function fEg. 

Imagine now that the frequencies correspond to 
differences of photon or phonon modes in a box of size L, 



and that 



then one finds 



27rri/L, n G Z, 

1 2rx 



Lr2 



G^^it)^~xr 



sinhr(<-L/2) 
sinhri/2 



(210) 
(211) 

(212) 



In the limit of infinite volume, linear response predicts 
that the retarded correlator falls off exponentially and the 
spectral function has the form of a Lorentziar^, 



-r\t\ 

ruj 



(213) 
(214) 



Non-linear effects in general dominate at large times with 
Git) ~ i-^/^ The cross-over time from the exponential to the 
power-law behavior is however parametrically of order A^^ in 
SU(iVc) gauge theories [lOl] . 
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Finite-volume spectral function 




Fig. 17. Comparison of spectral functions in finite and infinite 
volume, in a free (3+l)-dimensional massless scalar field theory 
at zero temperature. The infinite- volume spectral function is a 
step function at threshold, cu = |fc|. Two different volumes are 
shown. 



The lesson we draw from this example is that in infinite 
volume, the parameters of a transport peak could be read 
off the asymptotic i — ^ oo behavior of G^^{t). In finite vol- 
ume on the other hand, Eq. (j212L these parameters have 
to be read off for t <^ L, but large enough for contribu- 
tions from higher frequencies to have died off. The window 
of real-time in which the transport properties can be ex- 
tracted thus only exists for a sufficiently large volume. 
To conclude with this example, when P — > oo, p{uj) = 
the memory function becomes a delta function, f{t) — 
2^6{t), 7 = ^ and Eq. (|205p takes the form of the memory- 



free Langevin equation |100| . If one then attributes a ther- 
mal distribution to the occupation numbers of the pho- 
tonic energy levels, {a\ara) = (S„m"-B(w„), nB(w) = (e^"^- 
one recovers familiar results (see section |4?6)) such 

as 



1 
7rm 



doj 



(7/to)2 + u 



(215) 



where the dots represent a temperature-independent term. 
For 7/to — > 0, the first term reproduces the equipartition 
result ^UbT. 

To familiarize oneself with the infinite- volume limit 
further it may be useful to look at a specific example 
in 3-f 1 dimensions. Consider a free massless scalar field 
theory at T = and the correlator of the mass operator, 



GE{t,k)^ / dxe'"-'^ {(t,^{t,x)(j)\0)) 



(216) 



One easily finds the expressions for the Euclidean corre- 
lator 



GEit,k) 



-\k\t 



(47r)2i 

and the associated spectral function, 

1 



p{uj,k) 



(47r)2 



9iu~\k\). 



(217) 



(218) 



In a finite periodic box on the other hand, 

g-|p|t g-|fc-p|t 



GEit,k) 



-y 



2\p\ 2\k-p\ ' 



and the corresponding spectral function 



-P\ 



(219) 



(220) 



is, as expected, a collection of delta functions accompanied 
by a weight factor. The spectral functions are displayed in 
Fig. ([T7)) for illustration. The 'spikes' in the finite-volume 
spectral function have been drawn with a height propor- 
tional to the weight of the corresponding delta function. 
It is clear that the finite-volume spectral function does 
not converge towards the infinite spectral function point- 
by-point in w, but rather in the sense of distributions. In 
infinite volume the convolution of the spectral function 
with a Gaussian 'resolution function', 

F(wo,r)=/ dojp{u;,k = 0) ^- , (221) 

Jo V 27ri 

amounts to j^, for F ^ ojq. In finite volume, on the 
other hand, the corresponding integral amounts to (using 
the Poisson summation formula) 



(47r)2 



exp 



2 sin 



ujQ\m\L 
2 



a;oi|m| 



(222) 



r < Wo, r^L < cjo. 



Thus one sees that on a resolution scale F which is not too 
small compared to l/£, FL(ujQ,r) is a good approxima- 
tion to its infinite- volume counterpart. This point was al- 
ready made a long time ago in the study of non-relativistic 
systems [102] . The volume thus intrinsically limits the res- 
olution in u at which one can determine properties of the 
infinite- volume spectral function with good accuracy. The 
simple formula (j222p . displayed in Fig. ((TS]), also iUus- 
trates that the convergence of Fl to F is rather 'errati- 
cally' dependent on uq] for instance, if ujq is a multiple of 
2/i, all contributions of the type m — (0, 0, m) vanish. 



4 An Overview of Lattice Calculations 

In this section we give an overview of the technical aspects 
of lattice QCD relevant to the calculations of Euclidean 
correlators, and review the latest of these calculations; the 
discussion of their analysis in terms of the spectral func- 
tion is postponed till section ([S]). We begin by reminding 
the reader of the Matsubara formalism in section (|4.ip . 
The basics of lattice regularized QCD are given in sec- 
tion (14.21) . The question of the renormalization of corre- 
lation functions is briefly discussed in section (14. 3p . and 
their discretization effects are investigated perturbatively 



Harvey B. Meyer: Transport Properties of the Quark-Gluon Plasma 



25 



{Anf FlCmq, r) 



2 











rL = 2 

rL = 3 











30 40 50 60 70 80 90 100 

(Bp L 

Fig. 18. The function Fl as given by Eq. (|222[) . normahzed to 
converge to 1 in the hmit L — > cxd. It represents the convolu- 
tion of the finite-volume spectral function with a Gaussian of 
width r centered at luq. At fixed L, increasing the resolution 
in frequency comes at the cost of enhancing the finite-volume 
effect. 

in section (|4.4[) . We then review the lattice calculations of 
light-quark vector current correlators, HQET correlators 
and energy- momentum tensor correlators (sections 14.51 to 
14.71 respectively) ■ An important question for the field is 
how the accuracy of the calculations scale with the lat- 
tice spacing and volume, and this question is addressed 
in section (j4.8|) . Finally, we briefly mention the possibil- 
ity of probing the relaxation of conserved charge densi- 
ties through off-diagonal correlators, G^^ with A different 
from B. 



where 

Lsidtcl^it, x), <p{t, x)) - -C{idtc^{t, x),(t){t, x)). (227) 

For fermions, the representation of the partition func- 
tion using Grassmann variables leads to 

^fcrmion = /antipcriodic I^'/'] (228) 

eyig{- At J dxCE{dtiJj{t,x),'iJj{t,x))}. 

Without reproducing the proof in detail, we illustrate where 
the antiperiodic boundary condition comes from in a two- 
state system, {|0), |1)}. Coherent fermionic states are de- 
fined by 

|r;)=e-'"^'|0), = (0|e— '* . (229) 
where rj and 77* are Grassmann variables. The definition 

Jdr]* dr]r]r]* = 1 (230) 

furthermore implies the 'completeness relation' 

/dry*dr/e-"'"|77)(?7| = 1. (231) 

Some algebraic manipulations then lead to 

Tv{e-^"}^ (fc|e-^^|A:) (232) 

0<fe<2 

= Jd7fdT]e"^''^{-Tj\e~^"\7]), (233) 

where the minus sign in (— 77I translates into the antiperi- 
odic boundary condition. See for instance the appendix 
of Ref. [103) for the generalization to multiple degrees of 
freedom. 



4.1 Path-integral representation of the partition 
function 



4.2 Common actions in lattice QCD 



In this subsection we briefly sketch the origin of the imaginarj^or definiteness we write out the simplest action for lattice 



time path-integral formalism of thermal field theory. To 
arrive at the path integral formula for the partition func- 
tion (see Eq. I^Tj) of a scalar field theory, we evaluate the 
trace using the eigenvectors of the field operator 

${0,xM =^{x)\<j>). (223) 

Then the partition function can be written as 

Z = Tr{e-^"} ^JUJdM^)){Me-'''\M (224) 

We can then use the standard path integral formula for 
the matrix elements of U{t) (see for instance [53]). 

(MUma) = Ji^ldcb] (225) 
exp{-i/g dt' f dxC{dt(j>{t',x),(f>{t',x))}. 
Now setting t = — b ~ a and summing over a yields 

^boson = /peHodicM'/'] (226) 

exp{-/(f dt J dx CEidt(f>{t,x),4){t,x))}, 



QGD |104j . We restrict ourselves to a four-dimensional cu- 
bic lattice, and the lattice spacing is denoted by a. The 
unit vectors fl are aligned along the four directions of 
the lattice of highest symmetry. The dynamical variables 
are (a) the 'link' variables Uf^{x) S SU{3), which can be 
thought of as living on the link from 2: to x -|- afl, and (b) 
the Grassmann variables 'ip{x), which are Dirac spinors in 
the fundamental color representation and live on the lat- 
tice sites. The latter correspond to the quarks, while the 
'link variables' are classically related to the gauge poten- 
tial by L/^(x) = e'-^^f^^). 

The Wilson gauge action reads 



1 

7o 



ReTr{l 

P^iu{x) = Uf_,{x)U^{x + afl)Ui_,{x 
The Wilson fermion action reads 
Sf = a'^^-0(x) (Vp7p +h 



P^A^)}, (234) 
ai))-i/7^(a;)-i.(235) 



iflAj tPix), (236) 
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decade many improvements have been made to its earh- 
est implementation. Of particular relevance are the more 
faithful representations of the quark determinant through 
several pseudofermions (instead of just one), where the 
long wavelength modes are separated from short wave- 
length modes. The most widely used realizations of this 
idea are the mass preconditioning (MPHMC [117 ), the 
domain-decomposition (DDHMC [118, ) and the rational 
HMC (RHMC [nn]) algorithms. 

The calculation of correlation functions amounts to 
calculating quark propagators on background gauge fields, 
i.e. solving the equation D[U]ip = rj for a given spinor field 
r]{x) carrying Dirac and color indices. Much progress has 
been made in speeding up the solution of this equation, 
which results in an improved scaling of the cost towards 
the chiral limit [12T?1|121||122) . In the simplest implemen- 
tation, gluonic correlators are evaluated by computing the 
correlation on the generated ensemble of gauge field con- 
figurations. However, in the absence of fermions the lo- 
cality of the action can be exploited to sample the fields 
more efficiently (multi- level algorithms |1231I124) ). 

The 'quenched' approximation consists in neglecting 
the gauge-field dependence of the fermion determinant al- 
together. The fermions then only appear at the level of 
the correlation functions. As a consequence, the weight- 
ing factor of configurations e""^" is local, and a heat-bath 
update algorithm can be employed, resulting in an update 
speedup factor on the order of 100. See [125j for a review 
of the state-of-the-art in simulation algorithms. 

4.3 Renormalization of lattice correlation functions 



where the summation over fi is understood. Here = 
^(V^ -|- V*) is the symmetric covariant derivative and 
A = 9^9* the covariant Laplacian, 

V^^(.t) = - {U^,ix)4>{x + afi) ~ i;{x)) , (237) 
a 

V;^(a:) = i {^Pix) - U^{x)-'i^ix - afi)) . (238) 

In numerical implementations, the hopping parameter k, 
defined by 

^2am + 8, (239) 

is used to parametrize the bare quark mass, and the quark 
field is rescaled so as to give the 't/j{x)ip{x) term in Eq. (|236p 
a unit coefficient. The discretization errors of the Wilson 
action are 0(a), even at the classical level. Therefore it 
is useful to improve the action by adding dimension-five 
operators with coefficients tuned to remove the leading 
cutoff effects on the spectrum of the theory |105j . For- 
tunately, having made use of the equations of motion, 
a single operator needs to be considered, a Pauli term 
aCswi'{x)jaf_i^Ff^^{x)'ip{x) jl06| . The coefficient Csw = 1 + 
0((7o) has been determined non-perturbatively for a wide 
range of bare couplings in the iVf = fl07| , 2 [108] and 3 
theories |109l . 

The Wilson action (improved or not) preserves the 
vector flavor symmetry SU(iVf), while the axial symme- 
try is broken. The space-time symmetries are of course 
broken down to discrete groups on the lattice. The dis- 
crete symmetries P and C are preserved. In general, it 
is the exact gauge invariance of lattice QCD actions and 
the compactness of the gauge degrees of freedom make 
non-perturbative calculations possible. 

For a brief introduction to Kogut-Susskind (also called 
staggered) quarks [llOj . see [111) , section (3.F). An in- 
troduction to Domain- Wall pi^fTT^ and Overlap pH] 
fermions, which preserve chiral symmetry on the lattice, 
is given in [T1 I115] . 

Since the quark action is quadratic, the Grassmann 
variables can be integrated out exactly, yielding in the 
numerator the determinant of the Dirac operator D (or 
equivalently, oi Q = j^D = Q^). In the simplest form, the 
determinant is then represented as a Gaussian integral 
over a 'pseudofermion' field (j), which is a c-number field 
on which the Dirac operator acts in the same way as on a 
fermion field. The action reads 

5g[C/]+0t|Q|-^f0. (240) 

It is manifestly positive-definite for an even number of 
mass-degenerate flavors, while the absolute value of Q is 
taken for an odd number of flavors. The 'Boltzmann' fac- 
tor is a real, positive definite integrand in the path 
integral, and can therefore be interpreted as a probability 
distribution for the gauge fields. This property allows one 
to apply importance sampling methods to generate gauge 
field configuration which dominate the path integral. The 
best known algorithm to generate these configurations 
is the Hybrid Monte-Carlo algorithm [116] . In the past 



Most lattice calculations are restricted to on-shell corre- 
lators. An on-shell correlator avoids any contact terms in 
coordinate space, so that the equations of motion can be 
used to simplify its renormalization. In the case of the vec- 
tor and the axial-vector current only a finite multiplica- 
tive renormalization is then necessary. For instance, the 
renormalization factors ZA{go) and Zv{go) have been de- 
termined non-perturbatively in the Nf — 2 0(a)-improved 
Wilson theory [mi[T?7] . 

The renormalization of the lattice energy-momentum 
tensor is more complicated already in the continuum, due 
to the trace anomaly. In addition, the lattice regulator 
breaks continuous translation invariance, resulting in ad- 
ditional (finite) renormalization effects. In the continuum, 
a symmetric and gauge invariant form of the energy-momentum 
tensor is 

T^, = 0l^ + 0l^ + is^,{e^ + e'), (241) 

nS = 1 X pa pa _ pa pa 

_ <-> _ <-> _ <-> 



where D^—Dfj_ — D^. This form has finite matrix ele- 
ments between physical states (but divergences appear in 
offshell correlation functions). The term 9^ is the trace 
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anomaly [12811129) . On the lattice, there is a lot of free- 
dom as to how to discretize these expressions [130'. An 
important consideration is that the traceless part for 
fi = v belongs to a three-dimensional irreducible repre- 
sentation of the H(4) group |131) . while 9f_ii, ior n ^ v 
belongs to a six-dimensional representation and 6 belongs 
to the trivial representation. Because there are no other 
gauge invariant operators of dimension four in these rep- 
resentations, their renormalization is multiplicative. The 
renormalization of 9 is particularly simple when its dis- 
cretization is based on the action itself, in which case the 
renormalization factor is expressed in terms of the lattice 
beta function, dgQ^ /dloga. For {/j. = v), the renor- 
malization factor of 9)j,i, for n = v of the discretization 
based on plaquettes can be related to anisotropy coefh- 
cients P^IHg] . 

A number of quantities of interest are defined through 
frequency-space correlation functions, which are not on- 
shell quantities. These include the static susceptibility of 
the quark number density, 



Xs{k) = 




{n{-it,x)n{0)). (242) 



At fc = 0, this corresponds to the derivative of the free 
energy with respective to the chemical potential. Another 
example is the static susceptibility of the axial current, 
which was shown to be related to the finite-temperature 
dispersion relation of the 'pion' quasiparticle |134| . 

X5'5"'' = J dtdx{A-^{~it,x)A'a{0)), (243) 

A'^^ix) = hf.l5^i'- (244) 

Care has to be taken to renormalize these correlation func- 
tions. Using a conserved current greatly simplifies this 
task, but this requires an exact chiral symmetry on the 
lattice. Formulations preserving an exact chiral symmetry 
exist but are computationally expensive [135) . 

We illustrate in the case of the vector current how off- 
shell correlation functions on the lattice can be given the 
same structure as in the continuum |136] . By performing 
a non-anomalous change of variables in the path integral, 
one obtains the identity (SSO) — (SO) for a general oper- 
ator O. We now consider the case of Wilson fermions for 
simplicity. From the transformation 

S^{x) = uj{x)iIj{x), 5'4>{x) = -bj{x)-il>{x), (245) 

one obtains, for O — jy(0), 

Knt..{k) = 0, (fc^ - f sin(iafc^)), (246) 
7T^.(fc) = a4^e*-(j^(x)j.(0)) (247) 

X 



where the conserved vector current and the point-split 
scalar operator take the form 

j^{x) = \{^(x + ati){l+^^)U;;\x)i:{x) (248) 
- ■4j{x){1 - -i^)U^{x)'il){x + a^)). 

s,{x) - \{ij{x + av){l + -i,)U-Hx)^{x) (249) 
-I- ^{x){l - -iy)Uy{x)ip{x + av)). 

Eq. (j246p corresponds in configuration space to 

II^MOMla^b'.'W + a-Hf,J,,^s,{Q)) = Vx, (250) 

which shows that the second term takes care of removing a 
quadratically divergent contact term. No summation over 
V is implied in Eq. (PTTl) and 

4.4 Lattice perturbation theory 

The perturbative calculations of thermal spectral func- 
tions and Euclidean correlators are reviewed in section 
p.Sp . Here we describe the results of free field calcula- 
tions performed in lattice perturbation theory (see |138j 
for a general review of this subject). The motivation is to 
quantify the size of discretization errors for standard lat- 
tice actions and discretizations of the local operators in a 
controlled framework, namely leading order perturbation 
theory. Thus these calculations allow one to make an in- 
formed decision in choosing a discretization, and also allow 
one to estimate what lattice spacing it will be necessary to 
reach in order to achieve a certain accuracy. We will con- 
sider correlators in the mixed representation where they 
depend on time t and spatial momentum fc. This represen- 
tation has the advantage that the correlator is 'on-shell'. 
The approach to the continuum of on-shell correlators is 
considerably simpler than that of off-shell correlators be- 
cause the latter involve an integration over the contact 
region |105| . For a given lattice spacing a — l/{TNt), 
there will thus be a range of time separations t and spa- 
tial momenta fc for which the discretization errors is below 
a prescribed level, say 3%. 

As an example, it is interesting to quantify how badly 
the breaking of translation invariance on the lattice affects 
the conservation of energy and momentum as measured 
by a discretization of Too and Tofe- Indeed, the correlators 
Goo, 00(^7 fc = 0) and G03, 03(^1 fc = 0) are t independent in 
the continuum, but on the lattice they suffer of discretiza- 
tion errors which vanish as ~ for a fixed t. Fig. (IT9l) 
illustrates this point. It also shows that the discretization 
errors are numerically large, particularly in the energy cor- 
relator. A possible explanation for this large discretization 
effect is that the flattening results from a delicate cancel- 
lation of contributions which by power counting diverge 
as for small t. 

The resolution in the Euclidean time variable is ob- 
viously important to constrain the spectral function. But 
the computational cost of decreasing the lattice spacing 
grows very rapidly (see section (j4.8p . The situation may 
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Fig. 19. Top: the treelevel Goo.oo(a;o, 0) correlator at finite lat- 
tice spacing on the isotropic lattice in SU(3) gauge theory. The 
horizontal line is the continuum treelevel prediction, 47r^/15. 
Bottom: the treelevel Go3,o3(a^o, 0) correlator at finite lattice 
spacing on the isotropic lattice. The horizontal line is the con- 
tinuum trelevel prediction, 47r^/45 [137| . 



potentially be improved by choosing an anisotropic lattice, 
where the temporal lattice spacing is smaller by a factor ^ 
than the spatial lattice spacing. This strategy has been fol- 
lowed in certain thermodynamic studies |139) and in stud- 
ies of quarkonium in the deconfined phase of SU(3) gauge 



theory |140j . In perturbation theory it is observed that 
the discretization errors on thermodynamic potentials and 
on time-dependent correlators are more sensisitive to the 
value of the temporal lattice spacing than on the spatial 
one by a significant factor. At the same time the cost of 
producing independent configurations only scales as l/ur, 
whereas reducing the lattice spacing on the isotropic lat- 
tice would scale as 1/a*. This constitutes the main moti- 
vation to use an anisotropic lattice. An example is given 
in Fig. (PO)) . which displays the cutoff effect affecting the 
shear channel correlator G'i3.i3(i, fc = 0) in the SU(Afc) 
theory at maximum separation t = /3/2. Clearly, reducing 
the temporal lattice spacing overall reduces the discretiza- 
tion error, although it changes sign for a particular value 
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Fig. 20. The cutoff effects on the tensor correlator for xq = 
/3/4 and /3/2, corresponding to small and large symbols respec- 
tively [137| . The spatial lattice spacing is held fixed, the 
temporal lattice spacing a-r is varied between a quarter and 
one times a^- The D's refer to T12 and the (8)'s to ^(Tii —722). 
In addition, the two crosses at ^ = 1 indicate the reduction of 
the cutoff effect on the ^(Tn — r22) correlators when increasing 
/3/ao- from 8 to 16. 



of the anisotropy near 2. The figure suggests that reduc- 
ing the temporal lattice spacing by a factor 2 is almost as 
good as reducing the lattice spacing in each direction by 
the same amount. 

The results described above are obtained at tree level, 
and renormalization effects only set in at the next order. 
A drawback of the anisotropic action is that it requires a 
tuning of certain bare parameters in addition to the gauge 
coupling and the quark masses. In the pure gauge theory, 
the only such additional parameter is the bare anisotropy 
factor ^0 that enters the action. In effect it must be tuned 
in order to maintain a constant renormalized anisotropy ^ 
as the spatial lattice spacing is varied. The calculation of 
C(fo:5o) has been carried out using various methods for 
the Wilson action [14111142) . Including the quarks, there 
are additional dimension-four operators whose coefficients 
need to be tuned. Methods have been developed to achieve 
this in the context of hadron spectroscopy jl43j . In ad- 
dition the anisotropic lattice complicates the normaliza- 
tion of local operators used in calculating thermodynamic 
properties or correlation functions. A possible strategy is 
then to determine the renormalization factors of a par- 
ticular discretization based on the renormalization of the 
anisotropy ^; this discretization can thereafter be used to 
calibrate other discretization. 

The discretization effects in high-temperature meson 
correlators and spectral functions at nonzero-momentum 
and fermion mass were investigated by Aarts and Martinez- 
Resco in [SS] . In this case the spectral function can be an- 
alyzed straightforwardly at finite lattice spacing because 
the transfer matrix expression of the meson correlators is 
simple |103| . We give as an example the Euclidean cor- 
relator of a quark bilinear ipFip for (unimproved) Wilson 



Harvey B. Meyer: Transport Properties of the Quark-Gluon Plasma 



29 



fermions: 



GH{t,p) ^^Y.^r{S{t,k)rS{-t,p + k)r}{251) 



L3 



L3 



[ a^H^S,{t,k)Sl{t,p + k) 



(252) 



-a%^S,{t,k)Sl{t,p + k) 
-afs^{t,k)Si{t,p + k) 



where for < 7^ 



Si{t,k) = ^^""'"^l^^^^o^, cos\i{Ek{t - f )) , (253) 



2£kCosh{Ek/2T) 



S,{t,k) = - 



i sin ki 



e 2£fcCosh(£;fe/2T) 



sinh(£;fe(i- I)), 



Suit,k) 



■smh{Ek{t- f )) 



2£kCOsh{Ek/2T) 
and (r is the 'Wilson parameter' |104j . usuaUy set to unity) 



^Mk = r'^il - coski) 

3 

C/Cfc = ^ 7i sin fcj , 



cosh(£;fc/0 = 1 + 
£1 



2{l + Mk) ' 
{l+Mk)sm\i{EklO- 



(254) 

(255) 

(256) 
(257) 



Due to the UV cutoff, the spectral function on the lat- 
tice vanishes for frequencies larger than the energy of the 
largest momentum in the Brillouin zone, p = -^(1, 1, 1). 
This frequency is roughly 



L0„ 



2 , / 
-log 1 

a \ 



6- 



(258) 



In addition, the spectral function admits several cusps cor- 
responding to the edges of the Brillouin zone. For smaller 
frequencies, the finite lattice spacing affects predominantly 
the mismatch between the continuum and lattice light- 
cone, which can be substantial [66] . 

From a practical point of view, the conservative way 
to proceed is to extrapolate the Euclidean correlator to 
the continuum, see for instance |1441I145) . Then the ques- 
tion of the analytic continuation takes place entirely in the 
continuum |145) . This way of proceeding has the advan- 
tage that analytic studies of the spectral function in the 
continuum, where perturbative calculations can be pushed 
more easily to higher orders, can be used directly. 



4.5 The Electric Conductivity 



As reviewed in section p.2.ip . the electric conductivity is 
phenomenologically important in that it determines the 




Fig. 21. Vector current correlation functions calculated with 
Kogut-Susskind fermions on quenched lattices at zero spa- 
tial momentum. Top: T = 2Tc, with iVt = 8 (circles), 10 
(squares) and 12 (pentagons) [146) . Bottom: T ~ 0.6Tc with 
Nt — 24 |147l . On both plots the curves are the result of 
Bayesian fits. 



soft photon emission by the quark-gluon plasma, and it 
can be extracted from the vector current correlator via a 
Kubo formula. The early calculations of finite-temperature 
vector correlation functions were carried out with Kogut- 
Susskind fermions. The correlators from two such calcu- 
lations, by S. Gupta |146| (up to Nt = 12) and Aarts et 
al [I47] (up to Nt = 24), are displayed in Fig. As 
is seen most clearly from the bottom plot, there are ef- 
fectively two channels contributing, corresponding respec- 
tively to the vector current ip^iip and the quark bilinear 
'07o757i'0; the contribution of the latter comes with a -|- 
sign on odd time-slice and a - sign on the even ones |147) . 
Due to this effect, the analysis proceeds separately on the 
even and on the odd sites, effectively resulting in a loss of 
resolution in the time variable by a factor 2. For this rea- 
son most recent calculations have been based on Wilson- 
type fermions. 

A very recently published calculation using 0(a) im- 
proved Wilson fermions is displayed in Fig. (22), with 
temporal extents up to Nt 



48 [148]. The Euchdean cor- 
relator of the spatial components of the vector current 
at 1.45Tc has been normalized by its value in the non- 
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Fig. 22. Top: Vector current correlation function calculated 
with 0(a) improved Wilson fermions on quenched lattices at 
zero spatial momentum |145j and T ~ 1.45Tc. Bottom: from 
the same publication, the vector correlation function, Gii{T, T), 
normalized by the quark number susceptibility and the free 
vector correlation function calculated in the continuum (full 
symbols) and on lattices with temporal extent Nt (open sym- 
bols) . The results of five values of Euclidean time are displayed, 
tT = 0.1875, 0.25, 0.3125, 0.375 and 0.5 (top to bottom), on 
lattices with temporal extent Nt = 16, 24, 32 and 48. 
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Fig. 23. Vector charm current correlation functions in 
quenched QCD. The ratio to the reconstructed correlator 
from the confined phase is taken. Top: 0(a) improved Wil- 
son fermions on isotropic lattices at 1.5Tc [148) . The correla- 
tor in the denominator is reconstructed from 0.75Tc. In 'Vsub', 
the contribution of the transport peak has been subtracted. 
Bottom : anisotropic Fermilab lattice formulation (^ = 4) at 
several temperatures ( [140| . fig. from [149) ). The correlator in 
the denominator is reconstructed from a very low temperature. 
The enhancement over the reconstructed correlator is signifi- 
cant and is likely due to the contribution of a transport peak. 



interacting limit, as well as by the quark number suscep- 
tibility. The latter normalization factor has the benefit of 
cancelling out the normalization of the vector current. On 
the bottom panel, this data is extrapolated to the contin- 
uum. One notices that normalizing by the free correlator 
computed at the corresponding value of Nt reduces the 
slope of the extrapolation, without completely eliminating 
it. The figure thus shows that a continuum extrapolation 
is mandatory for a precision analysis of Euclidean correla- 
tors. The achieved accuracy is below 1%. The continuum 
extrapolated data provides a solid basis for a study of the 
corresponding spectral function, which we will discuss in 
section (|5.4p . 

The Euclidean correlator of the charm current has 
also been investigated extensively, with the primary mo- 
tivation to study the 'melting' of charmonium states in 
the quark-gluon plasma. This effect is expected on gen- 
eral grounds at sufficiently high temperature [150) . and 
is described in detail by resummed perturbative calcula- 



tions, see in particular |151[|152] . From the point of view 
of these studies, the presence of a transport peak of width 
^ r^/M, see Eq. ()137p . represents an undesirable contri- 
bution |153[|154] . which can approximately be subtracted 
by taking the difference of the correlator at two Euclidean- 
time points [T5511 148) . The functional form of the transport 
peak was calculated in [631 based on the Langevin equa- 
tion, and the result is reviewed in section p.3p . If the 
charm quark is heavy enough to be well described by the 
static approximation, then the parameters of the trans- 
port peak can be obtained from a correlator calculated in 
heavy-quark effective theory (see the next section). 

Figure (^5]) displays the ratio of the Euclidean cor- 
relator at l.STc to the reconstructed correlator taken at 
0.75Tc. The latter was obtained by first solving for p{uj, 0.75Tc), 
and then applying the definition ((58)) . The top panel presents 
data obtained on an isotropic lattice |148| , while the lower 
one uses an anisotropic lattice with aajar = 4 |140[ . These 
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Fig. 24. Thermal velocity of the heavy quarks as a function of 
temperature, as extracted from the area under the transport 
peak in the vector current correlator [155) . 



calculations are performed in the quenched approximation 
using 'clover' improved Wilson fermions. Similar studies 
had been carried out earlier [156] . An important result of 
these analyses is that once the transport peak contribution 
is removed from the high-temperature correlator, the sub- 
tracted correlator at 1.5Tc is indistinguishable from the 
correlator at 0.75Tc to the present accuracy (top panel of 
Fig. [231). 

At vanishing spatial momentum, the transport peak 
appearing in the density-density correlator is an exact 
delta function, and its weight is the static susceptibility 
Xs- The transport peak appearing in the current-current 
correlator on the other hand acquires a width of order 
T^/M. Its area is Xs{Vp) (see section [331) . Therefore the 
thermal velocity of the heavy quarks can be extracted 
from the ratio. The main task consists in separating the 
low-frequency from the high-frequency contribution in the 
current correlator. An interesting analysis was carried out 



in |155| assuming that for uj > 2M, the spectral density 
is unchanged from T = 0, which is consistent with the 
lattice data. Subtracting this high-frequency contribution 
from the Euclidean current correlator, what remains is an 
essentially i-independent function. Its ratio to the static 
susceptibility yields an estimate of the thermal velocity, 
displayed in Fig. ([24[) . These results are compared to a pre- 
diction based on kinetic theory. Equating expression (jl36l) 
with the susceptibility measured on the lattice defines an 
effective (temperature-dependent) quark mass [155 ; in- 
serting this effective mass into the equipartition prediction 
= leads to the black curve in Fig. ([24[) . Note that 
this curve is obtained without using the current correla- 
tor, and clearly the agreement with the data is excellent. 
This is evidence that the kinetic description works. The 
velocity appears to be solely a function of , rather than 
T and M separately. 

The transport peak has a width given by 77, the drag 
coefficient (see Eq. (|137l) ') which is of order 1/M . Because 
of the smoothness of the kernel relating the spectral den- 



9.4 



9.2 



C- 8.8 
8.6 



8.4 



8.2 



0.38 



A 



T=1.8T(, p=7.024 -^t- 
T=3.1Tc |3=7.483 
T=6.2Tc |3=8.090 
T=3.1T, 0(9") + 
Ak/T^=0.352 A=T X 



A 

+ 



0.4 



0.42 



0.44 
Tl 



0.46 



0.48 



0.5 



Fig. 25. Top: the Euclidean force-force correlator at f = 1/2T 
from simulations with the Wilson action at /3 = 7.483 and 
Nt = 8, 12, 16,22 |157j . The absolute normalization has been 
determined by matching the data at tT = | and T — Z.lTc 
with the NLO result of Burnier et al [65] . Bottom : the quan- 
tity n{t), defined in Eq. ([259)) . The denote the 0(g*) pre- 
diction at ?iTc [65] . The 'x' denote the function n{t) that 
results from adding the low-frequency contribution Apioj) = 
■irZiKtanh(tj/2T)0(/l— |w|) to the 0{g'^) spectral function, with 
parameters Ak and A given in the caption. 



sity to the Euclidean correlator, it is difficult to separately 
determine the width and height of the peak in the longitu- 
dinal current correlator from Euclidean data. This is cru- 
cial, since the height corresponds to the diffusion constant 
D. However, one can use heavy-quark effective theory to 
determine the momentum diffusion coefficient k, in terms 
of which D and 77 can be expressed as rj — k/ (2MT) and 
D = 2T^ / K in the heavy quark limit. This is the subject 
of the next section. 



4.6 Heavy Quark Diffusion 

In the current correlator of a heavy quark flavor, the tail 
of the transport peak contains the information on the mo- 
mentum diffusion coefficient k, see Eq. (jl39p . Here we re- 
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view a first calculation of the HQET correlator (11421) in the 
SU(3) gauge theory [157', from which k can be extracted 
from the Kubo-type formula (|143p . The classical interpre- 
tation of the correlator is that it measures the (Lorentz) 
force correlator along the worldline of a heavy quark. 

On the lattice the chromo-electric field is not protected 
from renormalization (unlike in the continuum). Since this 
normalizing factor has not yet been calculated, the next- 
to-leading order pcrturbative calculation of Burnier et al [51 
was used to normalize the lattice correlator. The condi- 
tion that the lattice correlator agree with the pcrturba- 
tive correlator at relatively short distances {tT = j) and 
T = S.lTc was used. The resulting values of the Euclidean 
correlator at the maximal separation of t = ^ are dis- 



100 



played in Fig. (j25p . The temperature was varied by chang- 
ing Nt from 8 to 22 at fixed = 7.483. The tempera- 
ture variation between 2Tc and 6Tc is on the order of 20%. 
The sign of the variation is the same as predicted by the 
perturbative expression: the magnitude of the force-force 
correlator in the gluon plasma increases as the tempera- 
ture is decreased. 

A quantity that measures the relative fall-off of the 
Euclidean correlator is the quantity f2{t) > defined by 



(259) 



GE{t-a/2) _ cosh [Q{t){/3/2-{t-a/2))] 
GE{t + a/2) ~ cosh [X?(t)(/3/2 - (t + a/2))] 

It has a continuum limit, in which J7tanh J7(/?/2 — t) = 
— ■^\ogGE{t)- One can interpret Q{t) as the location of 
a delta function in the spectral function which by itself 
reproduces the local fall-off of the Euclidean correlator. 
The function J7(i) is displayed in Fig. (|25l bottom panel) 
for three different temperatures, where the temperature is 
varied this time by varying the bare coupling at fixed 
Nt = 16. Because the statistical samples of the numerator 
and the denominator in Eq. (|259D are highly correlated, 
the numerical results for these ratios have uncertainties 
at the few-permille level. Figure [55] also displays the per- 
turbative prediction at 3.1Tc [5S]; it is based directly on 
Eq. (|259p rather than on the continuum version of this 
equation, to allow for a more direct comparison with the 
lattice data. The lattice J7(t) differs from the perturba- 
tive one only by a few percent, namely the latter falls off 
slightly more steeply. Since the difference is so small, it 
could partly be due to discretization effects. To reduce 
their infiuence, only the largest values of t were consid- 
ered. One may ask nonetheless, how large a difference in 
the transport coefficient k could this discrepancy possibly 
correspond to. An estimate is obtained by adding a low- 
frequency correction to the perturbative spectral function. 



Ap{u) 



1 



An tanh(a;/2T) e{A - \uj\). 



(260) 



Other functional forms than (j260p . such as a Breit-Wigner 
curve, would perhaps be more realistic, but would not 
change the conclusions in any significant way. Adding such 
a term to the spectral function has the effect of making 
the Euclidean correlator fiatter, and indeed, by adjusting 
An, one can obtain good agreement for the largest two 
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Fig. 26. Euclidean correlation function Ge{xo) of T\2 (top) 
and r,i^ (bottom) in the SU(3) gauge theory, at zero spatial 
momentum and several temperatures in the deconfined phase 
with Nt = 8 [14411159] . The solid curves represent the treelevel 
correlators. 



t values between the perturbative prediction modified by 
Eq. ([2601) and the lattice data. At T = S.lTc, for A = T 
and Ak/T^ = 0.352(38), agreement is obtained with the 
lattice data at the two largest t values. This represents a 
substantial enhancement of k over the leading-order per- 
turbative value. An equally good agreement is obtained if 
one chooses A = 2T, which leads to Ak/T^ w 0.204(22). 
While it is too early to draw phenomenological conclu- 
sions, this increase appears to be not quite sufficient to 
explain the experimentally observed elliptic ffow of heavy 
quarks [^DU^fTSS] . 

One of the applications of this type of calculation in the 
HQET framework is that it can provide the parameters 
of the transport peak in bottomium, and perhaps even 
charmonium correlators. We refer the reader to section 
(|3.3p for a review of this connection [64]. This could be 
instrumental in the analysis of the 'melting' of quarkonia 
states. 



4.7 Shear and Bulk Viscosities 

The calculation of the shear viscosity is of the greatest 
importance for heavy ion phenomenology. But for techni- 
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Fig. 27. The lowest two moments n = and 1 defined in 
Eq. 1551 normalized to their treelevel values 162 . For compar- 
ison, the same quantities computed in the strongly coupled 
A/" = 4 SYM theory are displayed [55] . 



cal reasons the associated spectral function is also proba- 
bly the hardest to constrain by using lattice simulations. 
Therefore results have only been obtained in the SU(3) 
gauge theory (see also the SU(2) results investigating the 
behavior of correlators near the second-order phase tran- 
sition [160] ). 

The Euclidean correlators in the shear and in the bulk 
channels are displayed in Fig. (j26p . They are given as func- 
tion of Euclidean time t and are evaluated at vanishing 
spatial momentum, fc = 0. The lattices used in these cal- 
culations are isotropic and have Nt = 8. The treelevel 
correlator is displayed as a solid black line. Clearly the de- 
parture from this curve is small. This is due to the large 
contribution from the ultraviolet frequencies, which re- 
sults in a behavior at short distances and dominates 
the Euclidean correlator for all t. The data was obtained 
with a two- level algorithm |1241I161) . which exploits the 
locality of the action and the operator T^i/, both of which 
are formed of color-traces of 1x1 plaquettes. The algorithm 
averages out the short- wavelength fluctuations of the two 
operators separately before evaluating their correlation. 
In this way, the 'statistics' accumulated for each of them 
multiplies with the statistics accumulated for the other. 
These stochastic sub-averages do not lead to a bias in the 
final expectation value |123| . As a consequence of using 
this algorithm, the relative statistical uncertainty on the 
data points is roughly independent of t. 

In order to see the departure from the free-field behav- 
ior, one has to inspect the correlator in greater detail. The 
moments (1551) are directly computable on the lattice. The 
first two moments of the shear correlator, which amount 
to the value and the curvature of the Euclidean corre- 
lator at t = J, are displayed in Fig. (|27p as a function 
of temperature. They were computed on A't = 8 lattices 
and are normalized to their treelevel values, which can 
be obtained straightforwardly by integrating the spectral 
functions (|16ip . While below T — 2Tc one observes a clear 
departure of the lowest moment from the non-interacting 
limit, a higher accuracy (both statistical and systematic) 
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Fig. 28. The Euclidean energy density correlator for fixed 
t = ^, as a function of the spatial momentum q |164| . In 
these calculations, the lattice size is 12 x 48'^ at the two lower 
temperatures and 10 x 20^ at the two higher temperatures. The 
lattice is anisotropic, aa/ar — 2. At l.OlTc the intercept is as 
high as 13.6(1.9). 



would be needed to see a statistically significant departure 
from the non-interacting limit at the higher temperatures. 
The next moment is overall closer to the non-interacting 
limit, as one would expect on the basis that it weights 
the higher frequencies more strongly than the lowest mo- 
ment. Before interpreting the departure of less than 10% 
from the free value, the discretization errors would need 
to be brought under a control comparable to the statisti- 
cal errors. Note, the strong statistical correlations between 
Euclidean data points result in a strong reduction in the 
statistical error for the second derivative. The figure also 
displays the ratio of the Af = 4 SYM moments, calculated 
at infinite and zero coupling ;56u80 . This confirms how 
numerically insensitive the Euclidean correlator is to in- 
teractions among the constituents of the plasma. It also 
shows that the size and magnitude of the deviations from 
the free approximation is the same in both theories. The 
numerical insensitivity of the Euclidean correlator to the 
transport properties was pointed out in |163) (weak cou- 
pling regime) and in [56 (strong coupling regime). 

The conservation of the EMT, d^T^^ = 0, implies in 
particular, for q = qe^ and with Euclidean sign conven- 
tions, 

P00,00(w, q) = P03,03{^, q) = g^P33,33(w, q), 

-t^^ Poi,oi(w, q) = pi343(w, q) . (261) 

In order to be more sensitive to the infrared part of the 
shear and sound channels, it is useful to combine the in- 
formation from all the Euclidean correlators that belong 
to these channels [ST] . Indeed the spectral function for the 
energy density correlator goes like at high frequencies 
instead of ijj^ . Similarly, the momentum density spectral 
function goes like uPq^ at high frequencies. 

The energy- density correlator at non-zero spatial mo- 
mentum q is particularly interesting. Its main contribution 
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Fig. 29. Computational cost scaling for the shear channel cor- 
relator [164] , 

comes from the sound peak, Eq. (|116p . The area under this 
peak is to leading order the specific heat, c^. At second 
order in (conformal) hydrodynamics, the position of the 
peak also depends on the relaxation time of the medium 
(Eq. 11231) . Figure (pSj) displays the Euclidean correlator at 
the fixed time i = f and as a function of q. The correlator 
falls off fairly slowly, by about a factor 2 between q = 
and |<7| — ttT. This is in qualitative agreement with the 
free field expectation. The main difference between the 
free field prediction, displayed as a black curve, and the 
lattice data is the value of the specific heat which gives 
the intercept in Fig. (^51) . In section (15. 6p . we describe an 
attempt to fit simultaneously several sound channel cor- 
relators as a function of both t and q. 

Very close to the phase transition, a much steeper fall 
as a function of q is observed. The phase transition is 
weakly first order in the SU(3) gauge theory. However the 
specific heat is very large just above Tc |165j . and this 
is confirmed by Fig. However for perturbations of 

wave- vector \q\ ~ ^ttT, the correlator is of the same typ- 
ical magnitude as at the higher temperatures. This sharp 
fall-off of the energy-density correlator between |q| = 
and \itT means that the area under the sound peak falls 
off rapidly in this wave- vector interval. This could be due 
to a corresponding fall-off in the static susceptibility of Tqo 
as a function of the wavelength of the perturbation. In any 
case, the Euclidean correlator is incompatible with the 
leading-order form of the spectral function (I116|) . Possi- 
bly, the fiuctuations of the order parameter (the Polyakov 
loop) have to be treated as slow modes in hydrodynam- 
ics (see for instance [166 ). We will return to the closely 
related question of the bulk channel near in section 
(EU). 



4.8 Computational Cost Scaling 

An important question for the field is how computation- 
ally expensive it is to improve our knowledge of the spec- 



tral functions. A first step is to look at the computational 
cost scaling of the Euclidean correlators. 

Let us consider the correlator of T12 in the SU(A^) 
gauge theory, out of which the shear viscosity can in prin- 
ciple be extracted. As long as we are in the deconfined 
phase, the correlator behaves as in a scale-invariant the- 
ory, up to small deviations. In other words, once the over- 
all dimensionality is scaled out by an appropriate power of 
the temperature T, the correlator is to a good approxima- 
tion a function of the variable tT. At the level of the lattice 
simulation, it is of course a function of Nt , Ns and the bare 
coupling (70 ■ The statement of scale- invariance translates 
into the lattice correlator being approximately indepen- 
dent of (70 ; as long as the system is in the high-temperature 
phase. The observable, normalized to be finite both in the 
continuum and infinite- volume limit, reads 

lim Iim^(/da;ri2(t,a;)/dyri2(0,y)\ . (262) 

L^oo a— >-0 h'^ \ I conn 

Using a straightforward algorithm to compute this observ- 
able, the root of the variance is approximately 

^L3(/dyri2(0,y)Ti2(0,y)) cx A^^^ (263) 

Li \ I conn 

Therefore the number of independent 'measurements' needed 
scales as N^N}^ . This is a very poor scaling indeed, and it 
means that reducing the lattice spacing by a factor two at 
a fixed temperature requires about three orders of magni- 
tude as much statistics. Note also that both the observable 
(|262p and the root of the variance are independent of the 
volume (in the large volume regime). 

The scaling is however dependent on the way the gauge 
fields are sampled. The fields can be sampled in two levels, 
and we denote by [. . . ] an average done at the lower level 
of the algorithm |123lll61j . Specifically, the estimator in a 
two-level algorithm reads 

( [/ dy ri2 (t , y)\ [/ do: T12 (0, a:)] ) . (264) 

Li \ I conn 

The root of its variance is approximately 

y2^L3/[/da;Ti2(0,a;)]'\ (265) 

Li \ I conn 

which scales in the same way as the observable. Therefore 
the number of measurements at the higher level of the al- 
gorithm is independent of A^t and N g. On the other hand, 
the number of measurements at the lower level of the algo- 
rithm must grow like in order to yield a good estimate 
of the sub-average [ / da; ri2(0, a;)] . Therefore the scaling 
of the number of measurements, all levels confounded, is 

Based on these power-counting estimates, the two-level 
algorithm yields a reduction of the number of measure- 
ments needed for a given statistical error by five powers 
of At. In both cases, the number of measurements must 
be multiplied by the volume, NfN^, to obtain the com- 
putational cost scaling. At the very least then, ignoring 
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any critical slowing down, the computational cost scales 
as NfN^ with the two-level algorithm. 

The prefactors of a scaling law are obviously impor- 
tant, too. Fig. (l29l) displays the actually performed num- 
ber of flops for a given statistical error in SU(3) gauge the- 
ory calculations |164) . These calculations were performed 
on anisotropic lattices ^ — 2, t is set equal to ^, and the 
parameter LT = Ng/Nt = 6 is kept constant along the 
curve. Therefore one expects a scaling with A^^^, which is 
consistent with the empirically observed variance of the 
observable, although the range of Nt is not long enough 
to exclude other scaling laws. 

A recent calculation of the isovector vector current 
correlator at l.5Tc in quenched QCD [145] shows that, 
at equal statistics, the error bars of the correlator at a 
fixed time separation only weakly depends on the lattice 
spacing. Since the configurations used were separated by 
a large number of update sweeps, one can assume that 
this reflects the behavior of the variance with the lattice 
spacing. It is clearly very different from the behavior in 
a flavor-singlet channel discussed above. This can be un- 
derstood on the basis of the fact that the variance can 
be interpreted in partially quenched QCD. The cutoff de- 
pendence of the variance is then dominated by the oper- 
ators of lowest dimension entering the operator-product 
expansion of the four-point function. The relevant opera- 
tor product is jiso(y)iiso(0), where jiso is the isovector cur- 
rent and jiso is a 'replica' thereof involving two unrelated 
flavors u and d. Because of the non-trivial flavor quantum 
numbers, no operator of dimension less than 6 can con- 
tribute to the OPE, and therefore the variance is UV-finite 
by power-counting. When the correlator is calculated with 
the 'point-to-air technique, neither the observable nor the 
variance have a volume dependence in leading order. The 
only source of increase in computational requirements as 
the lattice spacing is lowered is therefore the increase in 
the condition number of the Dirac operator, roughly oc Nt 
in the deconflned phase, and the obvious N^Nt increase 
in cost for handling a larger lattice. 



4.9 A variational method? 

One lesson of linear response theory is that the opera- 
tor B used to perturb the Hamiltonian need not be the 
same operator as the operator A, whose relaxation one 
wants to study. It is sufficient for it to have the same 
quantum numbers as A. One can even study the corre- 
lation of a whole set of operators B with A. Let us take 
A = n{t,k) = J da;e*''''^n(x) for illustration. Each of the 
correlators G^^{t) must decay with the same exponent 
exp— Z)fc^t, only the coefficient is given by the static sus- 
ceptibility = Gji^{ie), which depends on B: 

^AB t^oo ^AB^^-rt^ r{k)^Dk\ (266) 

The susceptibility xf^ expresses how well the operator 
B is able to induce an excitation of the particle number 
density with wavelength 27r/fc. If one were able to find 



an operator B such that G^^{t) is a pure decaying ex- 
ponential even at short times, then one could determine 
separately the static susceptibility xf^ and the decay rate 
r = Dk^. What functional form of the Euclidean correla- 
tor one should try to achieve by choosing an appropriate B 
operator is given by Eq. ([M)) . Since the exponential (|266p 
is the slowest decaying mode, the search for B amounts to 
finding the operator B whose Euclidean correlation with 
A decays as slowly as possible. This suggests a variational 
procedure to find the optimal B operator. Since this pro- 
cedure can be repeated for every fc, one can check that 
the decay rate F is quadratic in fc, thus confirming that 
one is in the regime described by hydrodynamics. We note 
that a variational approach to the calculation of charmo- 
nium spectral function has been proposed recently (in a 
different form) |167j . 

5 Euclidean Correlators and the Analytic 
Continuation Problem 

In this section, we discuss the problem of determining the 
spectral function p{lu) from Euclidean correlators. This 
question arises either when performing analytic calcula- 
tions in Euclidean space before continuing them to real 
frequencies; or when Monte-Carlo calculations are done 
in Euclidean space, and one wants to determine the main 
features of the spectral function. Unless explicitly indi- 
cated, we will restrict ourselves to the two-point function 
of a Hermitian operator. Therefore the spectral function 
is a real, odd function of frequency, and it is given by the 
imaginary part of the retarded correlator. 

As a historical note, this problem was first studied in 
the context of lattice QCD by Karsch and Wyld in the 
mid-eighties |168) . In non-relativistic systems, the prob- 
lem had already been formulated and different numeri- 
cal methods had been proposed to perform the analytic 
continuation [169 „l70|ll02j . A regularization of the inverse 
problem based on smoothness was introduced in [171) . The 
idea to impose a sum rule constraint was also used in |171) . 
In the early nineties, the maximum entropy method was 
introduced |172| and reviewed in |173j in the context of the 
many-body problem, as well as in |174) in the context of 
nuclear physics. The recent Ref. |175j compares different 
methods to extract the optical conductivity. There are still 
new methods that are being proposed [176] . see notably 
the stochastic method introduced in |177) . 

Let Gsi^^i) = be the frequency-space Euclidean 
correlator. It is initially calculated for the Matsubara fre- 
quencies cue = 27tT£. However, since it is analytic, it can 
be continued into the complex plane. In particular, based 
on Eq. the retarded correlator can be obtained by 

performing the substitution 

GRiio)^GE{uJe^-i[uJ + ie]), (267) 

and then using Eq. (|4T|) one arrives at the spectral func- 
tion, 

p{uj) = -ImGRiuj). (268) 

TT 
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Examples of spectral function calculations performed in 

]. The uniqueness 



this way can be found in |178l l 
of the solution is guaranteed by the Carlson theorem. This 
theorem of complex analysis implies that the only func- 
tion / that is analytic in the upper half plane, is bounded 
at infinity by e'*^^' with |c| < tt, and vanishes at z = in for 
n > is / = 0. The condition |c| < tt is obviously neces- 
sary, because sinhyrz fulfills the other conditions. The an- 
alytic continuation, as a mathematical problem, has been 
studied in detail in 11791. 



5.1 Formulation of the problem in configuration and in 
frequency space 

Here we will be mostly concerned with the numerical ver- 
sion of the problem. The first question one may want to 
address is in what form to formulate the problem. Both 
on the Euclidean side and on the Minkowski side, one has 
a choice of working in frequency space or in coordinate 
space. If one chooses frequency space on both sides, the 
naive procedure (based on Eq. I267p consists in fitting an 
analytic function to the retarded correlator Gr at the dis- 
crete Matsubara frequencies along the imaginary axis, and 
then evaluate the function for frequencies just above the 
real axis. One can also work entirely in coordinate space, 
in which case the kernel relating the Minkowski-time cor- 
relator and the Euclidean correlator is given in Eq. (|54p . 
One then makes an ansatz for the Minkowski-time corre- 
lator and fits it to the Euclidean data. The formulation 
of the problem that has been used almost exclusively is 
the formulation based on the coordinate-space Euclidean 
correlator GE{t) and the frequency-space Minkowski cor- 
relator, which are related by the relation (see Eq. ((5(I)) . 



GE{t) 



cosha;(/3/2 - t) 
sinhcj^/2 



(269) 



We are primarily interested in the low-frequency part 
of the spectral function. This leads to significant short- 
cuts. In the fully frequency-space based approach, the odd 
derivatives of the spectral function at the origin can be 
obtained from the derivatives of the retarded correlator 
along the imaginary axis, see Eq. (j53p . Since the trans- 
port properties can be extracted in this way via Kubo for- 
mulae, there is no need to literally evaluate the retarded 
correlator for an argument along the real axis. In the fully 
coordinate-space based approach, the low-frequency be- 
havior of the spectral function is encoded in the late- 
time exponential fall-off of the Minkowski-time correla- 
tor, which allows one to extract for instance a diffusion 
coefficient (G(t, k) ~ e~^^ *) without having to explicitly 
perform the passage to frequency space. Finally, in the 
formulation (|269|) . the transport properties are read off 
the spectral function using the Kubo formulae. 

Working in frequency space, one thus has to inter- 
polate between the Matsubara frequencies to obtain an 
estimate of the slope of the retarded correlator at the 
origin, yielding the transport coefficient. If one uses the 
coordinate-space Euclidean correlator as a starting point. 



one has to solve an integral equation numerically, either 
Eq. (I269P or Eq. ((54|) . In all three cases, the task repre- 



sents a numerically ill-posed problem. In the interpola- 
tion case, one may see this by noticing that in practice, 
only frequencies up to a certain magnitude are numerically 
available. In that situation, a polynomial which vanishes 
dX ijj = a ■ 2-kT for £ = 0, . . . , £max can always be added 
with an arbitrary coefficient to a given tentative solution 
of the problem without spoiling the agreement with the 
available data. The slope at the origin is then modified by 
an arbitrary amount. Of course, such a high-degree poly- 
nomial can be excluded if one takes into account the in- 
formation that the retarded correlator grows at most with 
a certain power (in an asymptotically free theory, or in 
theories satisfying non-renormalization theorems such as 
A/" = 4 SYM, this power is given by dimensional analysis). 



5.2 A numerically ill-posed problem 

To illustrate the ill-posed nature of the problem more 
quantitatively, let us consider the problem at zero tem- 
perature, when it reduces to the well-studied problem of 
inverting the Laplace transform. 



GE{t) - / A^e-^'p{uj). 
Jo 



(270) 



Bertero, Boccacci and Pike '180' proved the following prop- 
erty. The input 'data' Gf''''{t)^ GE{t) + SGsit) is nec- 
essarily imperfect, but is assumed to satisfy the quality 
criterion 



dt \6GE{t)\^ < 



(271) 



Let Sp be the function of frequency of which SGE{t) is the 
Laplace transform. Then the size of the set {6p}, in the 
sense of the L^(0,oo) norm, is infinite. Suppose however 
we introduce a limit on the 'oscillations' of 5p{lj) in the 
form 



dio uj\5p\uj)\'^ < 



(272) 



Let S be the set of 6p satisfying this bound. Then a lower 
bound on the L'^-size of S can be given for small e/E, 



sup5pesll^Pll 



> 



t:E 



2|log(e/£;)r 



(273) 



From here one clearly sees that (a) letting E be arbitrarily 
large results in an infinite L^-uncertainty on 6p and (b) 
for a given E, the L^-uncertainty only decreases logarith- 
mically in e. When Gf^'^it) comes from a Monte- Carlo 
simulation, improving the accuracy is thus exponentially 
computationally expensive in the desired L^-accuracy of 
the spectral function. 

While this fact may seem discouraging, the authors 



of |180| also prove that the situation improves drastically 
if the spectral function is known to satisfy certain analyt- 
icity properties. Suppose that p is analytic in the sector 
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we 



V|0| < a around the real axis, and that p is square- 



integrable along the direction 
estimate 



Then they provide the 



2a 



2a 



(274) 



In words, the uncertainty decreases much faster in this 
case. As a very favorable example, in the J\f = 4 SYM 
vector channel spectral function p69p . a — j. 



5.3 Linear methods 

We start by describing a class of linear methods to solve 
Eq. (|269p . A Fredholm equation of this type can triv- 
ially be modified by multiplying the kernel K{t,uj) = 

cash bJ {13 /2-t) 



[ihaj/3/2 



by a function of t times a function of w. 



GE{t) 



Auj p{uj)K{t, uj), 



(275) 



GE{t) = cj){t)GE{t), p{u)=p{u)/m{u), (276) 
k(t,cu) = (j){t)m{Lo)K{t,uj). (277) 



Since Eq. p75l) has the same form as Eq. (j269p . we will as- 
sume in the following that this 'preconditioning' has been 
done and that the kernel is regular everywhere. The twid- 
dle over Ge, P and K will be dropped. 

We assume that the input data in the inversion prob- 
lem is a set of N points of the Euclidean correlator, {Gi = 
GEiti)}fLi, along with their covariance matrix 



(278) 



For any linear method, the estimator of the spectral func- 
tion can be written in the form 



N 



(279) 



'fudged' version of the true spectral function. As the num- 
ber of points N increases, the resolution function becomes 
more sharply peaked, as illustrated below. 

As the simplest approach, one may choose a set of 
N linearly independent functions ui,{uj), and write jo as a 
linear combination of them, with coefficients to be deter- 
mined. 



N 

p{uj) = ^ atut{uj). 

£=1 



We also define 



Giit) 



duj K{t, uj) ui{uj) 



(282) 



(283) 



The coefficients ai in Eq. (|282p can be worked out ex- 
plicitly. The condition that p reproduces the Euclidean 
correlator amounts to 



N 

E 



Gi, 



(284) 



where we have set Gu = Gi{ti). 
ceed is to invert the matrix Ga 
decomposition, 

G = UwV^ 



A standard way to pro- 
using the singular-value 

(285) 

where U and V satisfy U^U = VV* = V*V ^ 1, and 
w — diag(ti;i, . . . , wn)- We assume the wi to be ordered, 
wi > W2 > ■ ■ ■ ■ One then obtains the solution of the linear 
problem. 



N N ^ 

qi{uj) = E E Vuiui{uj), 



£=1 t'. 



. A^,(286) 



Eq. (12861 12791) provide an estimator p of the spectral func- 
tion which reproduces exactly the input data points Gi 
when inserted into the integral equation (j275p . 



where the qi{uj) form a set of functions that depends on 
the specific method. Ignoring for now the statistical fluc- 
tuations affecting the values Gi, we can substitute for it 
expression (|275D into Eq. 



(I279D and obtain the relation 
between the estimator p{lo) and p(a;). 



/■CO 

p{uj) = / duj' 5{uj^uj')p{uj')^ 
Jo 

N 



(280) 
(281) 



Clearly, in view of Eq. (j280p the goal is to have the res- 
olution function d{uj,uj') as close as possible to a delta 
function (by a measure which is relevant to the problem). 
The finiteness of N results in a loss of 'sharpness' in re- 
lation (|280l) between the estimator and the genuine spec- 
tral function. The estimator p can thus be thought of as a 



5.3.1 Regularization 

So far we have ignored the finite accuracy with which the 
Euclidean correlator is known. The next question to ad- 
dress is therefore how accurately the Euclidean correlator 
GE{t) must be known in order to determine the coeffi- 
cients Of with good accuracy. If Su' is the covariance ma- 
trix of the Euclidean correlator points, then the covariance 
matrix of p{uj) and p(w') reads 



N 

a{uj,Uj')= ^ q^{uJ)Si^'q^'{LLl), 



(287) 



with qi{uj) given by Eq. (|286p . The problem that arises 
is that the W£ fall off very rapidly as the mode number i 
increases. In practice, the direct application of Eq. (|286l 



then typically leads to a wildly oscillating function p 
with 100% uncertainty. 
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It is the coefRcients of modes of higher ^ that are hard- 
est to pin down. Therefore the simplest cure to this insta- 
bility problem is to restrict i to run from 1 to M < N in 
Eq. (|282p . Since there are now fewer unknowns than in- 
put data, one performs an ordinary, linear fit to a lattice 
correlator, where the reads 

N 

xM = I] (G. - G,){S-%AG^' - G,0, (288) 

(289) 



Gi = duj K{ti,uj)p{uj). 
Jo 

However this straighforward way of proceeding is some- 
times unsatisfactory, because the transition (as M is in- 
creased) between having an unacceptably large value 
and having a poorly conditioned system (leading to large 
uncertainties in p) is quite rapid. Therefore it is prefer- 
able to have a continuous parameter which allows one to 
balance the stability of the procedure against a faithful 
description of the Euclidean data. This is achieved by min- 
imizing a function obtained by adding a regulating term 
to the (PMl) . niuhiplied by a 'Lagrange multiplier' A: 



:f^x' 



XR. 



(290) 



The minimization of with respect to the fit parameters 
ae can now be interpreted as the minimization of the 
under the condition that R take a given value. For A — 
l/jJL, it is also equivalent to the problem of minimizing 
R-\- fJbX^ , which means minimizing i? at a fixed value of the 
X^ . There is a significant amount of freedom in choosing 
R. 

If one has an a priori idea of the functional form of 
the spectral function, then R can be chosen to measure 
the departure of p from that functional form. Typically 
the expected functional form is smooth, so that the regu- 
lating term R tends to make the estimator p less wildly 
oscillating. This kind of regulator is used also in non- 
linear methods, such as the Maximum Entropy Method 
(see Eq. I316P . A typical example would be a 'kinetic en- 
ergy' term R = duj [-^Yg{uj), where g{uj) > is a 
weight function. 

A different kind of regulator directly aims at reducing 
the statistical variance of p, for instance 



R 



dco 



du' D{uj, uj')a{uj, Lo') 



(291) 



for some symmetric weighting function D{u,u'). If the 
goal is to minimize the variance of p at the origin, one 
chooses D{uj,uj') — S{ll!)6{uj'). This sort of regularization 
is employed e.g. in the Backus-Gilbert method (see [181) 
and Refs. therein), and has also been proposed in the con- 
text of the analytic continuation of Euclidean correlators 
in [T77]. 



space, one should first specify a scalar product, which is 
given by a measure. In t-space, the inverse covariance ma- 
trix {S~^)iii provides a natural measure, see for example 
expression (j288p . In frequency space, we define a scalar 
product via a symmetric kernel D{lo,lo') as in Eq. (|29ip . 
We will see that this kernel ultimately allows one to ad- 
just what part of the spectral function one would like to 
reduce the statistical uncertainty on, and also what shape 
one wants the resolution function to have. 

We define the symmetric frequency-space kernel 



N 



i/(w,w')- J2 {S'\^'K{t,,u;)Kit,,,u;'). 



(292) 



There are functions with positive eigenvalues that satisfy 
the following generalized eigenvalue problem. 



duj' H{uj,uj')ue{uj') (293) 



— wj I duj' D{ijj,uj')uf {uj') , 
Jo 

and their orthogonality and completeness properties take 
the form 



diJ Ui{u})D{uj,uj')ui^{uj') = 6u' , (294) 



J2 M^) / dw" D{uj', oj")ue{uj") = 5{u: - w')- (295) 
1=1 -^0 

We order the eigenfunctions so as to make the sequence 
LOi > W2 . . . monotonically decreasing. The motivation for 
choosing this basis is that the functions with largest 



contribute most to the correlator, as seen in Eq. (j298p be- 
low. Therefore it is natural to determine the components 
of ^{uj) = ^^^■^aiui{u}) in the linear subspace that the 
correlator is overall most sensitive too. 

With Gii = Gi{ti), one easily checks that 



N 



W, 



^w^} ^ Gu{S ^)ii'Gi'i' — Sec 



(296) 



In other words, if 5 = LsLl the Cholesky decomposition 
of 5, the N X M matrix 



^ -1 1 
Uu = y {L )ii'Gi'e — 
— We 



(297) 



is orthogonal, [/*[/ — 1. We now restrict ourselves to 
M = N. The solution to Eq. ((^M)) then yields a spec- 
tral function of the form Eq. (j279L with 



N ^ I 

(w) = ^{S~^)u' ^ — ■ Gi'iui{uj) 



(298) 



5.3.2 Choosing an appropriate basis of functions 

For concreteness we give a specific realization of the gen- 
eral procedure outlined above. Both in i-space and in uj 



The resolution function (j28ip is given here by 



6{U!, Uj') 



N 

E 

1=1 



Ui{uj) I du;"D{u;',u}")ui{uj"). (299) 
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Fig. 30. Non-vanishing eigenvalues wi corresponding to oper- 

for TV = 8 



ator H{lo,lo') (Eq. [305] with K{t,Lo) 
and Tti = I - ^, i = . . . Af - 1. 



If Siii is the covariance matrix of the Euclidean corre- 
lator points, then the covariance matrix of the coefficients 
ai is given by 

^ -<5«'. (300) 



Since the eigenvalues appear in the denominator of 
Eq. poop , the contribution of mode £ is guaranteed to be 
well determined (in terms of statistical uncertainty) if wi 
is not too small. Finally, one can convert the covariance of 
the coefficients ai into the covariance between two values 
of p, weighted by the kernel D, 

/>00 /"OO ^ -1 

duj duj' D{uj,uj')a{uj,uj') ^y^—- (^d) 
"'0 Jo 



To interpret Eq. (I299P and (|30ip . it is easiest to set 
D{uj,uj') — f{uj)S{uj — uj'). Then these formulae become 



N 



S{u;,uj') = J2M^)fi^>e(.^'), (302) 
e=i 

/ du;f{cj)a^{u:,u:) = J2—- (303) 

Thus one sees that if / is chosen to be peaked at some fre- 
quency the resolution function S{u!o,l^') will be more 
concentrated around that point. But this comes at the cost 
of the whole variance X^fci ^ being concentrated around 
luq. Suppose that the spectral function is expected to have 
a narrow peak at wq. To balance the two effects against ea- 
chother, one can maximize the figure-of- merit S/ -i/a^ + B, 
where S is the area under the peak, a = cr{uJo^ ^o) and B is 
the area under the 'smooth' part of the spectral function. 
The task is reminiscent of choosing a bin size in experi- 
mental particle physics when trying to extract a signal S 
from a background B. 
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Fig. 31. Resolution function 5(cj',aj) = '^f^-^Ui{ijj)ui{ijj') 
at the origin {uj' — 0), for different values of A'^. The 
itf(tj) are eigenfunctions of H{uj,uj') (Eg. 13051 with K{t,uj) — 

cosh uj{^ — t) ^ 



To regularize the problem, an option is to minimize 
+ Ai?, with R given by Eq. (|291l287p . but with 



N 



N 



q,{uj) = ^(5-1),,- ^z, • G,,iue{uj), (304) 

i' = l 1=1 

and the zi are treated as the fit parameters. 

5.3.3 Illustration and an example from the literature 

For illustration, a possible set of orthogonal functions is 
formed by the N non-zero modes of the kernel 



N 

E 



K{U,u)K{U,uj'). 



(305) 



Upon discretizing the frequency variable, Wj = jSuj, these 
modes can be obtained conveniently from the singular- 
value decomposition of the matrix Kij = \/SLijK{ti,ujj), 
if* = UwV^. 

The eigenvalues wi are displayed in Fig. ([50)) for the 
case m(a;) = tanh(a;/2T). Note that these eigenvalues are 
now solely a property of the kernel K{t, u). Unfortunately, 
they fall off very rapidly as i increases. This means that 
already the fourth or fifth normalized eigenfunction ui{u;) 
make only a small contribution to the Euclidean correlator 
(recall Eq. ([296])), in other words the correlator has little 
sensitivity to it. 

The resolution function centered around the origin is 
displayed in Fig. ([5T]) . It is clear that it becomes more 
strongly peaked as N increases, but that it takes quite 
a large N for it to become narrow on the scale of the 
temperature T. 
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OD/T 

Fig. 32. The spectral function p[ui,T) in the bulk channel 
of SU(3) gauge theory reconstructed by a linear method from 
Nt = 12 lattices at two different temperatures [144| . The bulk 
viscosity is estimated by C,/T^ = (tt/IS) x intercept. The oscil- 
lating curve is the (rescaled) resolution function 5(0, a; ). 



In [144] ■ the bulk channel was studied, more precisely 
the p = correlation function of 0. The spectral function 
was parametrized as 



with 



p{uj) 



m{uj) 



N 

m(uj)[l + ^ aiUi{uj)] 
e=i 



tanh(4a;) tanh^(c/?aj) 



(306) 



(307) 



The constant A is given by leading order perturbation 
theory, Eq. (11611 The sensitivity of the final result to the 
parameter c was probed by varying it by a factor 2 around 
the value j. The ui{uj) were the eigenmodes of the opera- 
tor (|305p . with K{t,uj) — m{u!) '^°^^"J'J and m{uj) given 



by Eq. p07p . Since N was only 4, the determination of the 
coefficients ai did not require invoking a regulating term. 
The resulting spectral function, as well as the resolution 
function associated with this reconstruction, are displayed 
in Fig. ([5^ . Due to the large UV tail of to(w), the reso- 
lution function 6(0, to) is rather broad in units of temper- 
ature. This clearly signals that only the gross features of 
p{i-u) are reliable. The actual spectral function may well 
contain more rapidly varying features which are washed 
out in this analysis, due to the paucity of the data. 



5.3.4 The method of Cuniberti, de Micheli and Viano 

The analytic structure of Minkowski and Euclidean corre- 
lators was studied thoroughly in P^79 . The basic quanti- 
ties defined there are the one-sided sum 



. t 

J+(t) 




(308) 



i>0 



Fig. 33. The complex planes, basic functions, and analytic 
structures (figure from [182] '). 



which is analytic for Imr < but has cuts in the upper 
half-plane; its discontinuity across the cut starting at the 
origin, 

J+it,-) = i[g+ie + it,-)-g+{-e + it,-)], (309) 

t > 0, e = 0+, which equals the retarded real-time corre- 
lator, TZ{t, •); as well as its Laplace transform 

J+(C,-)= rdte-<*J+{t,-), (310) 

which is analytic for Re^ > 0. For ( = w„, J+(C, •) re- 
duces to the Fourier components Q {ujn , • ) j and therefore 
constitutes the desired analytic continuation to a complex 
half-plane. The value of J^(C, •) along the axis Q = e — iuj, 
w S M, yields the Fourier transform of the retarded corre- 
lator, TZ{uj, •), whose imaginary part in turn equals the 
spectral function. The basic analytic structure is illus- 
trated in Fig. [221 

Cuniberti et al. |179| also proposed a practical pro- 
cedure to analytically continue the Euclidean correlator. 
It was recently reviewed and tested in a realistic exam- 
ple by Burnier et al. |182| . Here we follow the notation 
and treatment of the latter publication. From our point of 
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view the most important result of Cuniberti et al. is that 
the retarded correlator G{t) can be obtained as a series in 
Laguerre polynomials (normalized so that -^^(O) = 1), 



G{t) = exp [ - exp(-27rr<)] ^ at (26"^''"), (311) 
2(-l)' E —/-G'^e^'^ 2Fi(-^, n + 1; 1; 2)(312) 



n=0 



Here 2F1 is the hypergeometric function; Burnier et al. 
provide a recurrence relation to compute the combina- 
tion appearing in Eq. pi2p . From here the spectral func- 
tion can be obtained by Fourier transformation of G{t), 
Eq. p2p , exploiting if applicable the symmetry properties 
of G{t) Eq. (|28l29p . However, if the goal is to determine 
transport properties, they can also be read off directly 
from the i — >■ 00 behavior of G{t) (see e.g. Eq. 12131 and 
I214p . Finally we note that formula (|311l) is valid for Gsit) 
continuous at i = mod /3. In relativistic theories this 
property will in general only be fulfilled after a subtrac- 
tion of the hardest UV contributions. 

Since the required input into formula (j31ip is the Eu- 
clidean correlator in frequency-space, we may assume that 
the latter has been calculated on the lattice for the TV first 
Matsubara frequencies. The practical recipe that the re- 
sult pi2p suggests is to truncate the sum pill) at some 
^ = ^max |182j . Various criteria can be used to motivate the 
choice of £ma,x- In the original proposal, based on mathe- 
matical arguments the suggestion was to look for a plateau 



' as a function of imax- An alternative, more 
physical condition put forward in '182 is to require that 

G{t) |35], which is equivalent to X^fci = 0. We 
remark that through the dependence of ^max on the data, 
one introduces a non-linearity of the spectral function on 
the input data. 

In the case that the channel under consideration ad- 
mits a diffusion pole, the asymptotic behavior of G(t) is 
(disregarding non-linear effects that generically lead to 
power-law tails |101j ). 



Git) ^ xDk 



2 g-Dk^t 



(313) 



It is interesting to see how this behavior arises from a for- 
mula such as piip . At large times, e"^'^"^* will be very 
close to zero, and therefore the Laguerre polynomials are 
evaluated near the origin in this regime. Since the leading 
term, ^^a^, must vanish for the retarded correlator to 
vanish at infinity, the exponential fall-off must arises from 
the derivative of the Laguerre polynomials at the origin. 
It is thus clear that the expression with a finite number of 
terms will always fall off as e"^'^^* for sufficiently large t. 
However, this asymptotic behavior can of course be differ- 
ent when all the terms arc included. Specifically, we note 
the formula 



r(/3 + i)£(-i)^ 



j=0 



13 



U{x), (314) 



which for x — 2e~^'^"^* and (3 — Dk'^ /2ttT allows one to 
express the diffusion constant D directly in terms of the 
large-^ behavior of ai. 

Burnier et al. performed a practical test of the method 
on a UV-subtracted Euclidean correlator. With N equally 
spaced data points in Euclidean time, they tested the ac- 
curacy of the reconstructed spectral function, assuming a 
relative accuracy of 10^"^ or better on the Euclidean mock 
data. In this case, for the kind of test spectral function 
used, the accuracy achieved on p{uj) is better than 50%. 
Another lesson is that doubling N from 24 to 48 only 
moderately improves the result, and the accuracy must be 
drastically improved to benefit from the additional points. 
We expect the latter to be a general qualitative conclusion. 



5.4 The Maximal Entropy Method 

The Maximal Entropy Method (MEM) was first applied 
to the determination of spectral functions in lattice QCD 
in |183j . The method is described in detail in |184) (see also 
the older review [173] ) . Some technical improvements have 
been introduced more recently (' [18511186] ). It is notewor- 
thy that the method is also currently being used to study 
the spectral function in non-relativistic systems [1871. 

From a pragmatic point of view, this method is also 
based on the minimization of + ^^^^ where R acts as a 
regulator. But the difference with linear methods is that 
neither the nor R are any longer quadratic in p. Using 

a more standard notation the quantity that is maximized 

1 2 

is thus not simply e^s-**^ , but rather 



P cx cxp \ ^'2^ ^ '^^ 



s[9] = r 

Jo 



did 



p{uj) — miio) — p(i^) log 



p[uj) 



m{uj) 



(315) 



.(316) 



For Gaussian distributed data, minimizing the corre- 
sponds to the condition of maximum likelihood of the pa- 
rameters, given the observed data. The functional S is 
called the Shannon-Jaynes entropy, and measures the de- 
parture of the reconstructed p{u)) from the prior function 
771(0;). The relative weight given to the 'prior knowledge', 
encoded by 5, and the observed data, encoded by the x^j 
is fixed by the parameter a. We shall return to the ques- 
tion of how to set its value shortly. 

Similar to the linear methods, the singularity of the 
kernel should be removed by multiplying it by a function 
which is proportional to w at small frequencies, and com- 
pensating by dividing p{ijj) by this function 185,186 . 

The spectral function /o(w) is assumed to be positive. 
It is parametrized as follows. 



'p{uj) — m{uj) exp/(o;) 



(317) 



where m{uj) is the same 'prior' function as in Eq. pi6p . 
The exponent / is parametrized using a similar basis of 
functions as described in (15. 3p . 



M 



/(w) = ^aeue{u;) 



(318) 
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Fig. 34. The spectral function p{uj,T) in the vector channel 
of quenched QCD obtained by a MEM analysis [185j . The nor- 
malization of the spectral function differs by a factor 2n from 
ours; the conductivity a is 21 x the intercept. The prior used 
for p/u> was m{oj) oc {bNt + 



where are the M eigenfunctions of highest eigen- 

value of the eigenvalue problem 



((293)) . with M < N. 
In the recent literature, the basis was obtained by per- 
forming the singular-value decomposition of the matrix 
Kij = K{ti,u!j), K* = UwV*. This corresponds to choos- 
ing = J2tiK{U,uj)K{U,w') and D{uj,uj') - 
S(uj — Lj') in Eq. ()293D . The solution to the minimization 
of ~ ck'S' is then obtained iteratively, for example with 
the Marquardt-Levenberg algorithm |181) . 

A reasonable value of the parameter a is one for which 
the per degree of freedom is of order unity. More re- 
fined ways of choosing a, or more generally averaging over 
can be found in 



it, can be found in |184| . To calculate the error on the 
spectral function, a simple and robust approach is to use 
the jacknife method. Splitting the data into batches has 
the advantage of giving some information on the spread 
of jO, and has been found to yield somewhat more sta- 
ble results |187] . One may also quote the error not on a 
particular value of p{to), but rather on its average over 
a chosen frequency interval |184) . which has a stabilizing 
effect. An example of a spectral function in the vector 
channel of quenched QCD obtained from a MEM analysis 
is reported in Fig. (IMl) . See the recent reference [145) for 
further examples. 

The choice of non-quadratic regulator pi6l) has an 
information-theoretic justification. An argument for this 
form of S is given in [184j . The argument is based on the 
Bayesian idea that one has an a priori expectation for the 
outcome of p, which is encoded in the function miuj). The 
idea is then to construct the most probable spectral func- 
tion by maximizing the conditional probability P\p\D H], 
where D denotes the data and H the prior information. 

Generally speaking, non-linear regulators have an ad- 
vantage over linear regulators when the spectral density 



sis of functions, and therefore a modest increase of this 
linear combination can lead to a strong increase in p{uj) 
itself. For this reason, the maximum entropy method is 
probably the best currently available method to analyze 
quarkonium correlators, where rapid changes of the spec- 
tral function are expected around threshold. 

5.5 Frequency-space method 

Here we describe an approach to determining transport co- 
efficients proposed in 58j starting from the Euclidean cor- 
relator in frequency space. The idea is based on Eq. ([55)) . 
which relates derivatives of the retarded correlator in the 
imaginary frequency direction to derivatives of the spec- 
tral function along the real axis; and which relates 
the retarded correlator for frequencies along the positive 
imaginary axis to the Euclidean correlator. By the Carlson 
theorem, the exact knowledge of the Euclidean correlator 
at the Matsubara frequencies, where it can be calculated 
by Euclidean path integral methods, determines its an- 
alytic continuation to all frequencies, since physical cor- 
relation functions can grow at most with a power-law at 
large frequencies. At the point w it then coincides with 
the retarded correlator evaluated at iui, so that 



exhibits sharp peaks over a slowly varying 'background' [181 
This is quite easy to understand in the case of MEM, since 
it is the exponent in Eq. (|317p that is expanded in a ba- 



dGsi^^) 1^=0 dp 



duj 



duj ' 
d^p 



etc. 



(319) 
(320) 



doj^ duj^ ' 

This allows one, in principle, to determine the transport 
properties, since the right-hand side of Eq. (I319P is pre- 
cisely the expression appearing in the Kubo formulae. To 
illustrate this point, it is useful to note that the contribu- 
tion of a diffusion pole p = \ x—T-PrT^rrr^ to the Euclidean 
correlator takes the form 



XsDk' 



Dk^ 



(321) 



Thus one sees that a non-zero value of p{uj)/lj near the 
origin translates into a non-analyticity of the Euclidean 
correlator in oj^ around the origin. 

For a function known at uj^ — 2t:T£ = ujm ■ one can 
write a polynomial of degree n going through the n -I- 1 
first points {G^^}"^q. It has a slope at the origin given by 



G^(w = 0)eff 



(-1)^ G 



G 



(0) 



£\{n-ey. 



(322) 



When expressed in terms of dispersion integrals as in Eq. ([5B1 - 
57l) , the degree of convergence in the U V of this expression 
improves with the order of the polynomial: with some com- 
binatorics one can show that it is asymptotically p{uj) /oj^^^ 
for n even and p(uj) / oj^'^'^ for n odd. For instance, the 
slope at the origin of the n — 2 polynomial is given by 

1 



(1) 



2ujM 



-6a; 



M 



{3Gf - 4G 



(323) 
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where the hnear combination (j57p appears. Thus the slope 
of the quadratic polynomial is negative-definite, just as 
the actual slope of the retarded correlator -^Gii{iuj). One 
might be tempted to use it as a Euclidean estimator for the 
transport coefficient (for instance the diffusion coefficient 
D). However, if the spectral function admits a diffusion 
pole of width Dk^ <^ ujm-, expression (I323p amounts to 
-|Xs/wm, while actually -^GR{iLo) = It is clear 

why that is: a polynomial interpolating between the Eu- 
clidean points is only a good approximation to the re- 
tarded correlator if the latter does not contain a scale 
ojo Dk^ <C ujm over which it varies significantly. 

However, one knows a priori that the number density 
correlator must exhibit a diffusion pole. Therefore it is 
appropriate to 'help' the determination of the slope by in- 
cluding a term of the form (j32ip in a fit to the Euclidean 

data set {G^"*}. The remainder can be assumed to analytic 
in near the origin. We note that the idea of parametriz- 
ing the frequency-space Euclidean correlator by a rational 
function, and then obtaining the spectral function from its 
discontinuity across the real axis was one of the first pro- 
posed ^170j . Extracting a derivative along the imaginary 
axis should be comparatively more stable. 

No attempt has been made yet to use this method in 
practice, mainly because the calculation of off-shell corre- 
lators such as is technically more complicated than 
working in coordinate space. The attractiveness of this 
method is that one does not need to reconstruct explicitly 
the whole spectral function, and this is a special feature 
of the transport properties - as opposed to, for instance, 
determining the width of a quarkonium state. Recently, 
the method has been tested to determine the optical con- 
ductivity in a Hubbard model and found to work quite 
well in this case f |175) . see section [^TT]) . 

5.6 Combining Numerical Euclidean Data with 
Analytic Results 

In attacking the inverse problem of determining the spec- 
tral function, one can take a 'purist' approach, and apply 
a ready-to-use formula which is exact in the limit of per- 
fect Euclidean data, such as piip or p22p . At the other 
extreme, if one has control over the functional form of 
the spectral function, then a fit with a small number of 
free parameters is appropriate. However in practice and 
for the foreseeable future the first approach is bound to 
be of limited use due to imperfect data, while in QCD 
it seems unlikely that one will have good enough control 
over the functional form over the whole real axis of fre- 
quencies. Therefore it will remain useful to exploit prior 
knowledge about the spectral function, such as the para- 
metric location of a diffusion or sound pole, even when a 
'purist' approach is adopted. Here we wish to collect the 
kind of analytic information one may want to make use 
of. 

Knowing even partially the functional form of the spec- 
tral function is of course a great advantage. Information 
about the functional form is provided by 'straightforward' 



l.SSTc 2.32re free gluons A = oo SYM 

{■n+jO/s 0.20(3) 0.26(3) oo ^ 0.080 

2-KTTn 3.1(3) 3.2(3) oo 2 - log 2 « 1.31 

0.40(5) 0.51(5) 0.17 0.38 

Table 1. Sound channel spectral function at 2.3Tc. In the 
Table, lattice results are compared to free gluons [188| and to 
the strongly coupled SYM results [16ll59j . Stat, errors only are 
given. We expect C, to be negligible at these temperatures. 



perturbation theory at frequencies lo ^ T. Power cor- 
rections are predicted by the operator-product expansion. 
In QCD at zero chemical potential, the leading power 
corrections at finite temperature are suppressed by four 
powers of the frequency, because that is the dimension 
of the lowest-dimensional gauge-invariant operators (the 
mass operator t/;^ comes accompanied by one power of the 
quark mass, since it is a chirally non-invariant contribu- 
tion). At low frequencies and momenta on the other hand, 
hydrodynamics is expected to predict the pole structure 
of the correlators correctly, see section p.2p . The main 
unknown is the domain filling the gap between these two 
asymptotic regimes. This interval is precisely the region 
that the Euclidean correlator at the first few Matsubara 
frequencies is most sensitive too. 

Sum rules can be used as a constraint on the spectral 
function, if the RHS of the sum-rule is known. One exam- 
ple is the difference of vector and axial- vector current cor- 
relators (Ea. ll97)ll99p . Because the difference of the vector 
and axial- vector spectral functions falls off rapidly at large 
frequencies, the difference should be easier to parametrize 
and be more sensitive to the low-frequency region. The 
shear and bulk channels also obey sum rules once the 
T = correlators have been subtracted, Eq. (|196p and 
(|194p respectively. In all these cases a subtraction is nec- 
essary in order to make the spectral integral converge. 
This entails the loss of the positivity of the function one 
is trying to reconstruct. This is clearly a drawback, since 
the positivity property has a stabilizing effect. It has to 
be decided on a case-by-case basis whether it is prefer- 
able to deal with non-positivity or with a large ultraviolet 
contribution to the signal. This is the major added diffi- 
culty one faces in relativistic field theories as compared to 
non-relativistic systems. 

Exploiting the fact that different correlators at finite 
(w, k) are related by Ward identities can help to constrain 
the unknown spectral function. This observation is par- 
ticularly useful in the sound channel, where the energy- 
density correlator is most sensitive to the sound peak, 
while the longitudinal pressure correlator is much more 
sensitive to the UV tail of the sound-channel spectral func- 
tion. 

Thanks to the known tensor structure of the vector and 
energy- momentum tensor correlators (see section [^751) . one 
can exploit a richer set of Euclidean data if one is willing to 
analyze two channels simultaneously. For instance, letting 
6'q be the polar angle of qr, one can exploit also momenta 
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Fig. 35. The spectral function obtained from the fit 
((325)1 [164] . 



that are not aligned along a lattice axis via [164 



Jo smh 



(324) 



I Psnd (w, q) cos^ e'q + Psh(w, q) sin^ 0q|. 



Here psh and psnd are the shear and sound channel spectral 
functions, using the notation psnd('^,'Z) = — Po3,03('^, 'Z) 
and pshi^,q) = —poi,Qi{^,q)- On the anisotropic lattice, 
the correct normalization of all components of the energy- 
momentum tensor involves many independent normaliza- 
tion factors. Apart from reproducing the correct thermo- 
dynamic properties, they must obey the Ward identities. 
This has been used to determine some of the normalization 
constants in the SU(3) gauge theory |164j . In the latter 
reference, an attempt was made to fit the sound-channel 
correlators with the ansatz 



Psnd — Plow + Pmcd + Phigh, 



where 



tanh(a;/2T) 



(325) 



1 -I- (Tio;^ 



Pmcd(w,(J,T) 2 2, ^ A 

-oj q tanh ( — ) 



(cj2-c2(q)q2)2 + (/-^^g2)2 1-H0-2W2 



tanh(a;/2T) 
Phigii(^,g,T) 
tanh(a;/2T) 



2rV2 + (tj2_g2 



.... , 9 , iU-' , 

^2^2 tanh-(— ) 



2dA 



2r^l5(47r)2 ' 



M2)2 ' 
(326) 



and Cs{q'^) is given by Eq. (|123D . In the present specific ex- 
ample 7 parameters are fitted to a total of 48 data points, 
obtained at T = 2.3rc on a 16x48^ lattice with anisotropy 
o-a/o-T — 2. The parameters are Fg^ai, (72, tij, ^, cr, M . Fig- 
ure ([55]) displays the reconstructed sound spectral func- 
tion based on lattice data at 2.3Tc, and the results for the 
fit parameters of interest are given in Tab. (j5.6p . In this 
fit, the contribution of the (5-function to the T33 correlator 
at <7 = according to (I187P was not properly taken into 
account. This neglect leads to an overestimation of the 



sound attenuation length Fs 



iv+c 

e+p 



(although the error 



bars on the T33 correlator being significantly larger than 
on the To3 correlator, the effect may be small). 

The fit yields an acceptable correlated value per 
degree of freedom, x^/^^-O-i = 1.6. An interesting point is 
that no prominent features are needed to fit the Euclidean 
data apart from the hydrodynamics structure at low fre- 
quencies and the perturbative UV tails at high frequencies. 
In this sense the result is reminiscent of the strongly cou- 
pled SYM spectral functions displayed in Fig. p^ . Based 
on the perturbative calculation [BSj, at 2.3Tc the bulk vis- 
cosity is expected to be negligible compared to the shear 
viscosity. The values of the shear viscosity are numerically 
quite small, and the 'relaxation time' T77 is quite well con- 
strained by the fit. This comes from the sensitivity of the 
Euclidean correlator to the precise location of the sound 
pole. We note that the ratio ^ is strongly constrained by 
the requirement of causality |189] . namely 



cl 



(327) 



The values of the parameters given in Tab. (|5.6p are close 
to but consistent with this upper limit. 

A different strategy was adopted in [SS] to interpret 
the bulk channel Euclidean correlator using a maximum 



Harvey B. Meyer: Transport Properties of the Quark-Gluon Plasma 



45 



amount of analytic information. The temperatures stud- 
ied were 1.02, 1.24 and 1.65Tc in the SU(3) gauge the- 
ory. The first step is to switch to the Euchdean correla- 
tor of the operator Oi, oc (T^^ -I- (Sc^ — l)Too) defined 
in Eq. (|192p . This removes the delta function at the ori- 
gin from the infinite- volume spectral function of T^^ . Sec- 
ondly, the difference between the finite-temperature and 
zero-temperature spectral functions was taken, which cor- 
responds to Eq. ([501) in terms of the Euclidean correla- 
tors. While this removes the bulk of the UV contribution, 
which represents an undesirable background when deter- 
mining the low-frequency part of the spectral function, it 
also introduces the (delta-function) contributions of the 
stable glueballs with negative weight. A property charac- 
teristic of SU(7V) gauge theories is that the mass gap m 
is numerically large, as compared to Tc- In particular, the 
subtracted spectral function Ap{uj, T) = p{uj, T) — p{bj, 0) 
is still positive for ut < m. The lightest two scalar glue- 
ball (G) masses are known quite well, and also the ma- 
trix element fa = (vacjO^IG) for the lightest glueball 
has been calculated [190lll91j . Therefore is would be de- 
sirable to compensate for these contributions by forming 
Ap{uj,T)+fQS{uj—m)+fQ,S{uj—m'). This linear combina- 
tion is then still positive up to the two-glueball threshold 
2m, which is about 3GeV. However, due to the imperfect 
knowledge of the glueball matrix elements, the subtrac- 
tion of the one-particle contributions was not carried out 
explicitly. 

The asymptotic UV behavior of the difference of spec- 
tral functions is given by the operator-product expansion, 
Eq. psip . Numerically, the OPE was assumed to be re- 
liable starting at w « 4.8GeV. The formula (|18ip then 
predicts Ap to be very small beyond that frequency. This 
sets a boundary value for the spectral function beyond 
which it essentially vanishes. 

The bulk sum rule in the form Eq. (jl94p provides one 
global constraint on the subtracted spectral function in 
terms of thermodynamic quantities, 

^ r^Ap.{c.,T) = 3{l-3cl)~4'—^. (328) 
e + pJo uj e+p 

Secondly the value of the subtracted Euclidean correlator 
at maximum time separation i = § was used. 



A simple pai'ametrization of Ap(co,T)/(a) s) 



rp4 



[GE(t,T)-GT{t,T)[ 



(329) 



1 



e+PJa 



duj Ap^{uj,T) 
T sinh ^ 



Expression (I329p turns out to be positive while expression 
p28p is negative. Since both receive the same contribu- 
tion from low frequencies, but the higher frequencies con- 
tribute practically only to (Eq. I328P , one infers that the 
subtracted spectral function must be positive at small fre- 
quencies and negative at higher frequencies. This can be 
quantified as follows. The frequency axis was divided into 
three intervals, defined by the separation points wi « m 
and W2 — 4.8GeV. In each interval, the subtracted spectral 
function is assumed to to be a quadratic polynomial in ui, 




Fig. 36. A smooth parametrization of the subtracted bulk- 
channel spectral function in SU(3) gauge theory, Ap^iuj, T) / (w- 
s), compatible with the bulk sum rule and the lattice correlator 
at t = i3/2[58). Recall that the ratio C,/s is given by -I x the 
intercept of this function. 



and the function is required to be continuous and differ- 
entiable at the boundary points. Beyond UJ2, it is assumed 
to be given by the operator-product expansion prediction. 
Based on the two pieces of information p28p and p29p . 
the parametrization of the subtracted spectral function 
obtained in [58^ is displayed in Fig. 

This calculation thus provides evidence for the deple- 
tion of the spectral density in the region of those glueballs 
that are stable or lie slightly above the two-particle thresh- 
old, and secondly for the appearance of a positive spectral 
weight for w < m. The analysis proceeded by maximizing 
the use of available analytic information on the spectral 
function. 



5.7 Comparison of different methods 

It is of interest to know how well different inverse problem 
methods compare, even if the outcome of such a compar- 
ison depends to some extent on the character of the spec- 
tral function one is trying to reconstruct from a Euclidean 
correlator. Such a comparison was recently done by Gun- 
narsson, Haverkort and Sangiovanni |175] . They started 
from the spectral function 



ct(w) 
+ 



_ irpiuj) 



Wo 



W2 



l + (c./A)2 1 + ((^ - e)/r2)2 

^ (330) 



i + {{u + e)/r2rii + {oj/r,r' 

a superposition of Lorentzians weighted by an overal fac- 
tor. The values of the parameters were chosen to be 



2ttT 



= 0.71, 



= 2.9, 



27rT 



2ttT 

7.1, Wi = 0.3, 



^ = 9.5, (331) 
W2 = 0.2 . (332) 
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Fig. 37. Comparison of different metiiods to reconstruct a 
known spectral function from the frequency-space Euclidean 
correlator. See the text for a description of the methods. Study 
and figure from [175) . The temperature is 0.067eV. 



in this example. To put this in perspective, Moore and 
Saremi [192' obtained the spectral function in the shear 
channel of SU(3) gauge theory and found a transport peak 
of width roughly F 0.05g'^T. For an intermediate 't 
Hooft coupling value of Xh — W — 12 thought to be rele- 
vant at a temperature of a few |193| , this corresponds 
to 2^ =0.09 .. 0.12. Thus the width (pT|) should be 
considered to be quite large and for such a smooth spec- 
tral function (on the scale of 2ttT) one would expect the 
inverse problem to be a fairly controlled procedure. We 
remark that there are also channels where the QCD spec- 
tral functions do not have kinetic theory peaks, for in- 
stance the tensor channel of the energy-momentum tensor 
or the transverse channel of the vector current. So a sub- 
tracted spectral function in a strongly interacting gauge 
theory may in fact bear some resemblance with this form, 
displayed as a solid red curve in Fig. ([57]) . 

After obtaining the frequency-space Euclidean correla- 
tor via (j52p and giving them uncorrelated Gaussian noise 
of relative standard deviation ao, four different methods 
were used to reconstruct the spectral function. The four 
panels of Fig. ([57]) display these results. We briefly sum- 
marize what these four methods are, and refer the reader 
to the original article |175| for the details. 

1. The Fade approximation is one of the oldest methods 
used in this problem [169) . The procedure consists in 
fitting a rational function ansatz to the Euclidean data 
(i.e. at imaginary frequencies), and then evaluating the 
imaginary part of the function along the real axis. 

2. the second method is very similar to the method de- 
scribed in section (|5.3p . except that it is based on 
Euclidean-frequency data. The spectral function is ex- 
panded in a set of functions which are eigenfunctions 
of K^K, if K is the kernel of the transform ([5^ mod- 
ified by a weighting factor to take advantage of the 
finite (effective) support of the spectral function. The 
eigenvalues of K^K become small rapidly and a cutoff 
rtjy on their number is introduced. 

3. the sampling method: this is similar to the Maximum 
Entropy Method, where however instead of using a de- 
fault model the spectral function is given by an average 
weighted by P[a\GE]- 

4. the Maximum Entropy Method with a flat prior dis- 
tribution of the parameter a and a featureless prior 
function m{ijj) displayed in Fig. ([57b *). We note that 
the data is divided into batches and the result of MEM 
in each of them is averaged over. 

See Fig. (1571) for the outcome of this test. The sensitivity of 
the results to variations in the parameters was also inves- 
tigated by doubling Fi and by reducing ao by an order of 
magnitude. The authors conclude that the Fade approx- 
imation is somewhat less reliable than the three others, 
while no clear winner was found among the SVD, MEM, 
and MEM sampling methods. For the spectral function 
studied, all three provide a decent approximation in the 
range of parameters they investigated. Using several meth- 
ods has the advantage of giving some indication of which 
features are genuine and which are artefacts of a particular 
method. 
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In addition, the method of extrapolating the Euchdean 
frequency-space correlator to w = described in section 
(|5.5p provides a rather good estimate of the conductivity, 
(t(0), in this example. This is illustrated in the inset of 
figure ([57h ). It is a good method if one is only interested 
in the transport part of the spectral function. 



6 Summary and Outlook 

We started this review by presenting the motivation to 
study the transport coefficients of the quark-gluon plasma 
in heavy-ion collisions and in cosmology. We then sum- 
marized the theoretical framework which underlies the 
determination of spectral functions from Euclidean-space 
Monte-Carlo simulations. We reviewed the analytic meth- 
ods that have been used to study the spectral functions 
in various kinematic regimes, such as the hydrodynamic 
expansion and kinetic theory at low frequencies, and per- 
turbation theory and the operator-product expansion at 
high frequencies. The recent lattice QCD calculations of 
Euclidean correlators in the vector, shear and sound chan- 
nels were reviewed in section 

In section ([S]), we summarized the main methods used 
to infer information about spectral functions from Eu- 
clidean correlation functions. Whatever conclusion one draws 
about the spectral function must be shown to be a prop- 
erty of QCD as opposed to a combination of QCD and 
model information. In most cases that were described in 
this review, the functional form of the spectral function 
is known both far in the ultraviolet lo ^ T (from pertur- 
bation theory) and deep in the infrared, w <C (from 
hydrodynamics), where tr is the longest relaxation time 
in the system. However practically no solid information is 
available for the frequencies between these two extremes 
for temperatures relevant to heavy-ion collisions. In addi- 
tion the onset of the validity of perturbation theory and 
of hydrodynamics is not known a priori, adding to the 
uncertainty of the ansatz one should use. 

We highlighted the differences between linear and non- 
linear methods to solve Eq. (I269p . The main advantages 
of the former is that it allows one to predict straightfor- 
wardly the maximum resolution one may possibly achieve 
in frequency space for a given kind of input data, before 
the data is available, and secondly that the solution is 
found directly by inverting a matrix rather than by an it- 
erative procedure. An attractive feature of the maximum 
entropy method (MEM) is that it is based on a generaliza- 
tion of the maximum likelihood principle to the situation 
where one has prior knowledge on the spectral function. 
In addition, it conveniently enforces the positivity of the 
spectral function by writing its estimator as the exponen- 
tial of a real function. The functional form of the ansatz 
means that one can more readily describe a spectral func- 
tion that contains both sharp 'peaks' and smooth 'back- 
grounds', but the difficulty of determining both the width 
and height of narrow peaks remains. 

In view of the inherent difhculty of any analytic con- 
tinuation, it is also worth trying to develop methods that 



circumvent this step, such as finite- volume and variational 
techniques. We have only briefly described some aspects 
of these research directions in Sections p.8p and (|4.9p . 
For instance, one should study what signature a diffusion 
process occurring in a finite volume leaves on current cor- 
relators. 

Nonetheless, what have we learnt about QCD spec- 
tral functions from the corresponding Euclidean correla- 
tors? Fit ansatze motivated by hydrodynamics at low fre- 
quencies and by perturbation theory at high frequencies 
are able to describe the Euclidean correlators, see for in- 
stance Fig. ([55)) where the estimated relaxation time was 
Tn — 3.1 /2nT at 1.58Tc. Recently, the isovector vector 
current correlator (extrapolated to the continuum) was 
successfully fitted by an ansatz composed of the free-field 
contribution, plus a Breit-Wigner peak centered at the 
origin. The relaxation time corresponding to the width of 
the peak is well constrained by the fit and about 5.6/27rr 
at 1.45Tc. These values certainly suggest that the relax- 
ation times are not long on the thermal time scale, as ki- 
netic theory parametrically predicts (l/g'^T) . On the other 
hand, they are not as small as their value in the Af ~ 4 
SYM in the strong coupling limit, where Tn ~ 1.3/27rT. 
The fact that the relaxation time is longer for the vector 
current is natural in perturbation theory, since the color 
charge carried by the quarks is smaller than the charge 
carried by the gluons 0. While the ordering of scales char- 
acteristic of a weakly coupled plasma is not respected at 
temperatures relevant to heavy-ion collisions (we knew 
this already from static quantities), some remnant of the 
kinetic theory transport peak would thus appear to be 
present, at least in the vector channel (see Fig.[Ml)- These 
observations should be regarded as preliminary; one's con- 
fidence in them would increase if tuT was seen to become 
progressively larger as the temperature is raised. By the 
methods used, it cannot be excluded at present that the 
longest relaxation times in the system are in reality much 
longer than these values, and that the spectral function 
has a more complicated structure than the fit ansatze al- 
low. The quoted results merely provide the 'simplest' ex- 
planation for the lattice data. 

One of the goals of this review is to indicate the prospects 
for decisive progress in the field using lattice QCD meth- 
ods. This exercise splits up into two separate questions. 
The first is, what accuracy is necessary on the Euclidean 
correlator in order to impact our knowledge of the spec- 
tral functions. And secondly, how much computing time is 
necessary to obtain Euclidean correlators with that preci- 
sion. We can illustrate this set of questions with the case 
of the force-force correlator that determines how fast a 
heavy quark diffuses in the quark-gluon plasma (see sec- 
tion ITH]). Addressing the first question, we want to de- 
termine the form of the spectral function at small fre- 
quencies, and it is important that one is really interested 

Comparing the widths of the transport peaks in Fig. ((S]) 
and (|10p . the ratio of relaxation times is approximately given 
by the ratio of the Casimir operators in their irreducible repre- 
sentations, Cf/Ca, which amounts to 4/9 for the SU(3) color 
group. 
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in a subtracted spectral function where the ultraviolet 
~ tail has been removed. An important aspect here 
is that even in the weak-coupling limit, no narrow fea- 
tures are predicted to occur in the spectral function. The 
analysis of |182| . where the narrowest feature of the (re- 
summed perturbative) spectral function has width « ttT, 
then shows that one needs a relative error at least as small 
as 10~^ on the UV-subtracted Euclidean correlator and 
Nt = 24 to reconstruct semi-quantitatively the spectral 
function. In this estimate, no physically motivated ansatz 
is made, it is a truncation of an exact mathematical for- 
mula. 

What amount of computing time is then necessary to 
determine the Euclidean correlator to this high degree of 
accuracy? Recently, the force-force correlator was calcu- 
lated in the SU(3) gauge theory on lattices with Nt be- 
tween ranging from 8 to 22 |157) . An Nt = 16 simulation 
around 6Tc cost about 2Gflop years (see Fig. [55]). The cost 
increases rapidly with Nt : to achieve the same accuracy at 
Nt = 22 as at Nt — 16, an increase in CPU-time by a fac- 
tor 20-25 is necessary. From this data, we estimate that it 
would take about half a month on a machine that sustains 
lOOTflops to achieve the relative accuracy of 10"'^ on the 
shape of the UV-subtracted correlator at Nt = 22. In ad- 
dition to the statistical errors, it is important to control 
the systematic errors at the same level. Thus one needs 
to check for finite- volume effects, and (probably more im- 
portantly in this particular case) for discretization errors. 
This requires varying the lattice spacing within a signifi- 
cant range. We note however that the overall normaliza- 
tion of the correlator does not need to be known to better 
than a few percent: it is the 'shape' of the correlator that 
must be controlled to very high accuracy in order to con- 
strain the spectral function. 

The example above shows that when detailed analytic 
information is available, and that relatively efficient nu- 
merical techniques can be used, the problem of determin- 
ing the spectral function and in particular the transport 
coefficient say at the 30% level can be tackled with cur- 
rently existing computing resources, albeit at the high end. 
It appears likely that the calculations will continue to be 
performed in the pure Yang-Mills theory for some time, 
since they are much faster and the crucial dynamics in 
any case appears to be associated with the gluons. Thus 
with a dedicated effort it is likely that in the next five to 
ten years, decisive progress will be made in this direction. 
Fortunately this is also the time scale on which heavy-ion 
collision experiments are approved for, in particular the 
CBM experiment at FAIR. The only assumption (backed 
up by resummed perturbation theory) made in this analy- 
sis on the spectral function is that it does not vary strongly 
on a frequency intervals of about ttT. Even before such 
a model-independent reconstruction of the spectral func- 
tion becomes possible, one may be able to draw interest- 
ing lessons already when the Monte-Carlo data quality is 
lower, for instance by making weak-coupling or strong- 
coupling (AdS/CFT) inspired fits. 

The isovector vector channel for light quarks is another 
case where we expect significant progress to be made, be- 



cause the cost does not increase as fast with Nt, and the 
spectral function increases more mildly, as w^, at large fre- 
quencies. On the other hand, the spectral function has a 
narrow transport peak at small frequencies in the weak 
coupling regime (parametrically, the width is g'^T and 
numerically, OAT to 0.6T for intermediate coupling val- 
ues (192) ). and determining separately the width and height 
of the transport peak is extremely difficult in that regime. 
The difficulty may be overcome with extremely accurate 
data, and possibly exploiting small but non- vanishing mo- 
menta by using twisted boundary conditions [194,195j. 

Given the observation of elliptic flow at RHIC 
[7] and most recently at the LHC ^ the shear viscosity 
is undoubtedly the most important transport coefficient 
for heavy-ion phenomenology. The shear and sound chan- 
nels are presumably the hardest to tackle by Euclidean 
methods, because they both exhibit a parametrically nar- 
row transport peak at weak coupling, they have a very 
hard UV contribution (w'* in the spectral function), and 
the variance is strongly UV-divergent. The multi-level al- 
gorithm [16V only partly alleviates the latter problem by 
reducing the power with which the variance diverges. So 
far these major difficulties have only partly been overcome 
by fitting a particular ansatz (motivated by hydrodynam- 
ics and perturbation theory) simultaneously to several cor- 
relators related by Ward identities (see section [S^ . 

The methods discussed here may be of interest to physi- 
cists in other fields where quantum Monte-Carlo simula- 
tions are used. In particular we have in mind cold Fermi 
gases, with which very controlled experiments can be made. 
When tuned to 'unitarity', the atoms have a divergent 
scattering length in the s-channel, and the system is strongly 
interacting. The shear viscosity to entropy density ratio 
has been estimated to be below unity, thus making this 
system an excellent fluid (see the review [5]). On the the- 
ory side, the spectral functions of a Fermi gas have re- 
cently been studied in kinetic theory [196) . An important 
technical simplification over the relativistic case discussed 
in this review is that the shear spectral function goes 
to zero at high frequencies |1971I198) . which implies that 
the Euclidean correlator is far more sensitive to the low- 
frequency part of the spectral function relevant to trans- 
port properties. 
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